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FOREWORD 


This  report  is  one  of  two  computer  program  user's  manuals  prepared  by 
Aerotherm  Division  of  Acurex  Corporation  under  USAF  Contract  F04611-69-C-0081. 
Included  herein  is  Volume  I  of  the  manual  for  Version  3  of  the  Aerotherm  Axi- 
Symmetric  Transient  Heating  and  Material  Ablation  (ASTHMA3)  computer  code.  This 
volume  describes  the  problems  solved  by  the  code  and  presents  an  input  (card 
format)  user's  guide  and  sample  problems.  The  code  was  originally  developed 
under  USAF  Contract  F04611-67-C-0047,  and  upgraded  under  the  subject  contract. 
The  work  was  administered  under  the  direction  of  the  Air  Force  Rocket  Propul¬ 
sion  Laboratory  with  Mr.  Robert  J.  Schoner  as  Project  Officer. 

Mr.  John  W.  Schaefer  was  Program  Manager  and  Mr.  Mitchell  R.  Wool  was 
Program  Engineer.  The  ASTHMA  code  upgrading  was  performed  by  Dr.  Carl  B.  Moyer 
and  Mr.  Kurt  E.  Suchsland. 

This  technical  report  has  been  reviewed  and  is  approved. 

A.  D.  Brown,  Jr.,  Lt.  Col.,  USAF 
Chief,  Technology  Division 


ABSTRACT 


This  document  presents  a  user's  manual  for  the  Aerothcrm  Axi-Symmetric 
Transient  Heating  and  Material  Ablation  Program,  Version  3  (ASTHMA3) ,  including 
a  general  description  of  program  capabilities  and  solution  procedures  as  well  as 
a  detailed  set  of  input  instructions. 

ASTHMA3  is  a  transient  heat  conduction  program  for  two-dimensional,  axi- 
symmetric  bodies.  Multiple  non-charring,  anisotropic  materials  may  be  studied. 
The  surface  boundary  condition  has  three  options,  including  an  Unusually  general 
thermochemical  erosion  or  ablation  condition  as  well  as  simplified  radiation 
and  specified  temperature  options. 
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GLOSSARY 


Some  of  the 
are  collected  here 

back  wall 

box 

column 

heated  surface 

mesh 

nodal  box 

nodal  center,  nodal 

node 

null  node 

row 

side  wall 


terms  used  in  the  text  have  very  particular  meanings.  These 
for  ready  reference. 


6ne  of  the  four  sides  of  the  nodal  mesh 
layout,  located  by  convention  at  the  bottoms 
of  the  columns;  may  be  convectively  and 
radiatively  cooled  but  does  net  ablate 

see  "nodal  box" 

refers  to  one  direction  of  the  nodal  mesh, 
heated  surface  is  at  the  top  of  the  columns 

the  one  surface  of  the  four  sides  of  the 
nodal  mesh  layout  which  is  exposed  to  the 
hyperthermal ,  chemically  reactive  environ¬ 
ment,  located  by  convention  at  the  top  of 
the  columns  in  the  mesh;  opposite  side  of 
mesh  network  is  the  back  wall 

assemblage  of  quadrilateral  nodal  boxes 
arranged  in  rows  and  columns,  each  row  and 
each  column  having  the  same  number  of  nodal 
boxes;  encompasses  all  the  material  of  in¬ 
terest  and  divides  it  into  nodes  for  finite 
difference  heat  conduction  analysis 

one  quadrilateral  zone  in  the  mesh 

point  a  point  within  the  nodal  box  at  which  all 

the  material  in  the  box  is  presumed  to  be 
lumped  for  the  finite  difference  calculation, 
not  necessarily  in  the  center  of  the  node 

used  for  nodal  box,  nodal  center,  nodal  point 

a  nodal  box  in  the  mesh  which  contains  no 
material,  used  for  bookkeeping  operations 
only,  not  involved  in  finite  difference 
solution 

refers  to  one  direction  in  nodal  mesh, 

"parallel"  to  heated  surface 

any  side  of  nodal  box  other  than  heated  surface; 
has  same  boundary  conditions  as  "back  wall" 


Preceding  page  blank 
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GLOSSARY  (Concluded) 


surface 


surface 


•iode  the  top  node  in  each  column  (excepting  null 

nodes) ,  that  node  next  to  (adjacent  to,  at) 
the  heated  surface;  note  that  the  surface 
nodal  box  is  next  to  the  heated  surface  and 
one  of  its  four  sides  is  exposed  to  the 
hyperthermal  environment,  the  nodal  point  or 
nodal  center  however  is  somewhere  within 
the  nodal  box;  there  is  a  special  point 
on  the  surface  called  the  surface  point 
which  must  not  be  confused  with  the  surface 
nodal  "center" 

point  the  center  of  the  heated  surface  side  of  the 

top  nodal  box  in  each  column  (not  counting 
null  nodes) ;  this  point  is  distinct  from 
the  nodal  point  for  that  nodal  box;  the  tem¬ 
perature  of  this  point  is  determined  by  the 
surface  energy  balance  in  Options  1  and  3 
or  by  assignment  in  Option  2 


LIST  OF  SYMBOLS 


A 


a,b,c 


B' 


b 


c 


CH 


c 

c 

6 


area  of  side  of  nodal  box 


nodal  box  side  and  center  line  lengtha,  Figure  8, 
Equation  (14) 

defined  as  (B '  -/,m  /p  u  CM)  thermochemical 
ablation  parameter  r£ 


defined  as  mc/peueCM,  total  ablation  parameter 

defined  as  /p„u  Cu 

JLd  ro  e  e  M 


defined  as  ™g/PeueCM 


same  as  B' 
see  a,b,c 

nodal  capacity.  Equation  (16) 

Stanton  number  for  heat  transfer 
(corrected  for  "blowing",  if  necessary) 

Stanton  number  for  heat  transfer  not  corrected 
for  blowing 

number  of  k  atoms  in  molecule  i 
Stanton  number  for  mass  transfer 
specific  heat 
see  a,b,c 

constant  defined  by  Equation  (39) 
binary  diffusion  coefficient 


ft2 

in 


Btu/°F 


Btu/lb°F 

in 

ft2/sec 
f  cVt>=c 
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LIST  OF  SYMBOLS  (continued) 


radiation  view  factor 

empirical  factors  appearing  in  Equation  (39) 
denotes  general  functional  relationship 
recovery  enthalpy 
enthalpy 

back  wall  convective  coefficient 

enthalpy  of  ablating  material  at  wall  temperature 

enthalpy  of  s  adjacent  to  the  wall 

total  number  of  identifiable  species 

number  of  iteration  cycles  in  problem 

chemical  species  indices 

total  number  of  elements  in  system 

mass  fraction  of  species  i 

mass  fraction  of  element  k  (regardless  of 
molecular  configuration) 

thermal  conductivity 

forward  rate  constant  for  kir.etic&lly 
controlled  reaction 

path  lengths  in  a  nodal  box,  Figure  6 


Btu/lb 

Btu/lb 

Btu/ft2sec  °R 
Btu/lb 
Btu/lb 


Btu/ft'-sec°R 

various 

ft 
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LIST  OF  SYMBOLS  (continued) 


£  total  path  length.  Figure  13 

Vi  system  molecular  weight  J}i*i 

molecular  weight  of  species  i 

m  node  corner  and  center  row  number 

m  mass  flow  rate  per  unit  area  from  the 

surface,  thermochemical  effects  only 


ft 

lb/lb  mole 
lb/lb  mole 

lb/ft2 sec 


m  total  mass  flow  rate  of  "char"  or  main 

ablating  material  per  unit  surface  area,  all  lb/ft2sec 

effects  (thermochemical  plus  condensed  phase 
mechanical  removal) 

m  mass  flow  rate  of  pyrolysis  gas  our  a  unit  lb/ft2sec 

g  area  of  surface 


flow  rate  of  condensed  phase  £  mechanically  lb/ft2 sec 

removed  from  surface 


N 


total  number  of  nodes 


n 

P 

Pi 

qchem 

qcond 

qdiff 


node  corner  and  center  column  number 
total  pressure 

partial  pressure  of  species  i 

chemical  energy  flux  term  defined  by  Equa¬ 
tion  (42) 

rate  of  energy  conduction  into  solid  material 
at  surface 


rate  of  energy  input  to  solid  surface  by 
diffusional  processes  in  the  boundary  la^er 


qrad  in 


rate  of  energy  input  to  the  surface  by  radia¬ 
tion  from  the  boundary  layer  or  from  outside 
the  boundary  layer,  same  as  q  ^ 


lb/ft2 

lb/ft2 

Btu/ft2sec 

~3tu/ft2sec 

Btu/ft2sec 

Btu/ft2sec 
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LIST  OF  SYMBOLS  (continued) 


grad  out 

^sen 

R 

R* 


rate  of  energy  radiated  away  from  surface 


4  <-eVB«r  -  »e  > 
w 

thermal  resistance,  see  Equations  >17) , (18) 


contact  thermal  resistance 

radius 

seeAS 


Btu/ft2sec 

Btu/ff"sec 

sec°F/Bfcu 

ft2sec°F/Btu 

in 


computing  time 


hrs 


temperature 


w 


U' 


Xi 


wall  (surface)  temperature,  general  term 

for  T  _ 
s ,  n 


thermal  conductance.  Equation  (20) 


modified  U,  Equation  (21) 


velocity  of  gases  at  edge  of  boundary  layer 


nodal  volume 


gas  velocity  (see  pv) 


mole  fraction  of  species  i 


axial  coovdinate 


Btu/sec°R 
Btu/sec0  R 
ft/sec 
ft3 

ft/sec 


ft 


zl 


diffusion  driving  force,  see  Equations  (36) 
and  (37) 
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LIST  OF  SYMBOLS  (continued) 


GREEK 


heated  surface  absorbtivity  (taken  equal  to  sw) , 
thermal  diffusivity 


Y 

A 


constant,  empirically  chosen  =  2/3 


denotes  change 


AS 

A0 


change  in  surface  point  location  during  A9 


time  step  in  finite  difference  solution 


emissivity  or  emittance 


input  multiplicative  safety  factor  in  time 
step  calculation.  Equation  (23) 


time 


constant  in  Equation  (41) 


dimensionless  factor  defined  by  Equation  (38) 


(pv) 


w 


density 

®c  -jy 


in 


P  u 
Ke  e 


m^ss  flow  at  boundary  layer  edge 


peueCH 


heat  transfer  convective  film  coefficient 


peueCM 


mass  transfer  convective  film  coefficient 


Stefan-Boltzmann  constant 


<P 


parameter  defined  by  Equation  (41) 


ft 

sec 


sec 

lb/ ft3 
lb/.  2sec 
lb/ft 2  sec 
lb/ft 2 sec 
?b/ft2.-5ec 


Btu/ft2sec°Rl* 


LIST  CP  SYMBOLS  (continued) 


SUBSCRIPTS 


A,B,C,B  indices  for  path  lengths  L  (Pig.  6) ,  nodal  side 

areas  A  (Fig.  7)  and  thermal  resistances  R  (Fig.  12) 


abs  denotes  absorbed  energy  rate  not  proportional  to  surface 

area;  Figure  15 


B  see  A,B,C,D 

bw  denotes  back  wall  side 


C 

CL 

c 

cond 

D 

diff 

e 

g 

H 

i/j 

k 

i 


see  A,B,C,D 

denotes  nodal  column  center  line 

denotes  "char"  or  ablating  material;  see  m  ; 
denotes  nodal  mesh  corner  coordinate  v 

denotes  heat  conduction;  Figure  15 

see  A,B,C,D 

denotes  diffusive  energy  fluxes.  Figure  15 

denotes  boundary  layer  outer  edge,  or  boundary  layer 
edge  gas  (environment) 

denotes  pyrolysis  jas 

see  C„ 

H 

denotes  any  identifiable  species;  atom,  ion,  molecule 
denotes  element 

index  of  condensed  phase  species  mechanically  removed 
from  surface 


M 


see  C., 

M 


m 

N 

n 

P 

r 

rad  in 

rad  out 

res 

s 

w 

r 

1,2 

0,0' 


' (prime) 
’ (prime) 


LIST  0?  SYMBOLS  (concluded) 

SUBSCRIPTS  (concluded) 

node  comer  and  center  row  number  index 

denotes  node  center 

node  corner  and  center  column  number  index 

denotes  virgin  plastic 

see  m 
r£ 

denotes  radiation  flux  to  surface;  Figure  15 

denotes  flux  radiated  from  surface;  Figure  15 

denotes  "reservoir"  communicating  with  back  wall  and 
side  walls 

denotes  heated  surface 

■v 

denotes  wall,  i.e.,  heated  surface 
see  C 

H0 

denote  "earlier"  and  "later";  also  see  Figure  13 
at  times  0  and  0',  respectively 
SUPERSCRIPTS 

sea  B',  3^,  B£,  B^,  B^,  U' 
denotes  at  new  time  8'  =  0  +  A0 
enthalpy  evaluation  temperature 
SPECIAL  SYMBOLS 

see  Z|,  R* 
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SECTION  1 


This  report  is  Volume  I  of  a  two  volume  user’s  description  of  the  Axi- 
Symmetric  Transient  Heating  and  Material  Ablation  Program  (generally  denoted 
ASTHMA) ,  Version  3  (ASTHMA3) .  This  volume  presents  an  analysis  technique  and 
describes  the  associated  computer  program  (ASTHMA3)  for  the  prediction  of  the 
thermal  response  and  ablation  of  two-dimensional  axi-symmetric  non-charring 
materials.  Section  2  of  this  volume  (Volume  I)  presents  some  of  the  underlying 
analysis  foundation  upon  which  the  computer  program  is  built.  Since  much  of  the 
pertinent  analysis  work  has  previously  been  published  elsewhere,  these  particular 
parts  of  the  analysis  are  treated  only  briefly  here  and  the  reader  is  referenced 
to  the  previous  reports  for  more  information.  Section  3  of  Volume  I  is  a 
complete,  but  condensed  and  succinct,  user’s  manual  giving  explicit  instructions 
for  the  preparation  of  input  to  the  program  and  the  interpretation  of  program 
output.  The  in-depth  temperature  prediction  is  of  the  familiar  explicit  finite- 
difference  type.  It  allows  a  completely  general  finite  difference  mesh  layout 
relative  to  the  physical  r-z  axes,  and  accounts  for  anisotropic  heat  conduction 
effects.  The  heated  surface  boundary  condition  is  an  unusally  general  thermo¬ 
chemical  type.  It  accounts  for  two  specific  kinetically  controlled  carbon 
oxidation  reactions,  one  specific  kinetically  controlled  carbon  reduction 
reaction,  one  kinetically  controlled  water  gas  shift  reaction,  and  any  number 
of  gas-phase  equilibrium  reactions  for  any  combination  of  ablating  materials 
and  environments.  Chief  applications  for  the  computer  program  are  rocket 
nozzles  and  entry  -'ehicles.  The  following  sketch  serves  to  clarify  various 
physical  aspects  of  the  ablation  problem  treated.  The  table  below  cites  the 

Physics  Computation 


boundary 
layer  edge 


Calculate  radiation  flux  to  wall 
Determination  of  boundary  layer 
edge  state 

Account  for  boundary  layer  trans¬ 
port  of  mass  and  energy  to  and 
from  surface 

Determine  chemical  state  at  sur¬ 
face  and  perform  energy  balance 
coupling  to  in-depth  solution 

In-depth  thermal  response 
calculation 


aspects  and  indicates  which  descriptions  of  the  analysis  are  to  be  found  in 
the  present  report  and  which  are  to  be  found  elsewhere.  Note  that  the 
ASTHMA3  Pregram  described  in  this  report  is  usually  used  in  conjunction  with 


SUMMARY  OF  COMPUTATIONAL  TASKS 


I 

Task 

1 

Program 

Where 

Described 

1. 

Determine  radiation  flux 
from  free  stream  to  wall 

Not  computed,  but  accounted 
for  by  ASTHMA 3 .  User  must 
provide  as  input. 

2. 

Determine  boundary  layer 
edge  state  'chemical) 

GASKET,  etc. 

Refs.  1-5 

3. 

Account  for  boundary  layer 
transport  of 

a.  mass 

GASKET,  etc. 

Refs.  1-5 

b.  energy 

ASTHMA,  but  user  must  pro¬ 
vide  convective  transfer 
coefficients  as  input 

4. 

Surface  thermochemical 
state 

GASKET/  etc. 

Refs.  1-5 

5. 

In-depth  temperature 

ASTHMA3 

Below 

another  program  which  handles  the  chemical  state  computations  required  in  the 
heated  surface  boundary  condition  definition.  Any  one  of  several  thermochemistry 
computer  programs  is  suitable  for  this  purpose,  including  EST  (Ref.  1) ,  EST2 
(Ref.  2),  EST3  (Ref.  3),  ACE  (Ref.  4),  and  GASKET  (Ref.  5).  These  codes  all 
differ  to  various  degrees  in  details  of  the  thermochemical  solution.  The 
GASKET  version  is  the  newest  and  contains  a  simplified  chemical  kinetics  treat¬ 
ment  of  carbon  oxidation  particularly  pertinent  to  rocket  nozzles.  Consequently, 
this  manual,  when  dealing  with  communication  between  ASTHMA3  and  the  surface 
state  solution,  will  refer  only  to  GASKET  as  the  partner  program.  Neither 
GASKET  nor  the  other  chemistry  codes  is  described  in  this  manual.  References 
1-5  provide  appropriate  descriptions  and  user's  manuals. 

A  complete  ablation  problem  may  be  solved  with  the  GASKET  and  ASTHMA 
programs  in  conjunction  in  the  following  steps: 

(1)  Lay  out  nodal  geometry,  provide  this  and  material  properties 
data  to  ASTHMA3 

(2)  Determine  radiation  flux  history  and  convective  heat  transfer 
coefficient  (peuecH)  history;  provide  to  ASTHMA3 


1-2 


(3)  Determine  chemical  nature  of  environment  and  generate  surface  thermo¬ 
chemical  ablation  tables  with  GASKET  for  input  to  ASTHMA3 

(4)  Run  ASTHMA3,  obtain  temperature  and  surface  recession  histories. 

The  second  volume  of  this  user's  description  presents  listings,  flow 
charts,  and  lists  of  FORTRAN  variable  names  for  the  program. 
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SECTION  2 


DESCRIPTION  OF  ANALYSIS  FOUNDATIONS  FOR  THE 
AXI- SYMMETRIC  TRANSIENT  HEATING  AND 
MATERIAL  ABLATION  PROGRAM,  VERSION  3  (ASTHMA3) 

This  section  describes  much  of  the  analysis  work  upon  which  the  ASTHMA3 
program  is  built.  Section  2.1  below  defines  the  problem  treated,  and  Section 
2.2  gives  some  description  of  the  in-depth  solution  procedure.  Section  2.3 
briefly  describes  the  heated  surface  boundary  condition  treatment. 

2.1  PROBLEM  TREATED  AND  GENERAL  METHOD  OF  SOLUTION 

2.1.1  Problem  Description 

The  basic  problem  is  to  predict  the  surface  and  in-depth  temperature 
history  and  surface  location  history  of  a  two-dimensional,  non-charring,  aniso¬ 
tropic  (but  orthotropic)  insulating  material  exposed  to  a  chemically  reactive 
hyperthermal  environment.  The  chief  practical  examples  are  rocket  nozzles  and 
nose  tips;  the  general  problem  is  depicted  schematically  in  Figure  1.  Below  the 
heated  surface,  the  material  response  is  characterized  by  anisotropic  heat 

_  _ Axis  of  Symmetry _ _ 

heated,  ablating  surface 


Figure  1.  Sketch  of  Ablating  System  Considered 
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conduction  with  temperature  dependent  thermal  conductivity  and  specific  heat. 

At  the  heated  surface  there  is  a  complex  boundary  condition  involving  thermal 
convection  from  the  hot  gases  adjacent  to  the  surface,  diffusive  mass  trans¬ 
port,  and  chemical  corrosion.  The  surface  chemical  reactions  may  be  in  equili¬ 
brium  or  may  be  kinetically  controlled  in  certain  respects  (e.g.,  the  oxida¬ 
tion  of  carbon  surfaces) . 

2.1.2  Summary  of  Solution  Procedure 

The  solution  procedure  chosen  for  this  problem  is  complex  and  is  de¬ 
scribed  in  more  detail  in  Sections  2.2  and  2.3  below.  As  a  quick  summary,  the 
in-depth  procedure  is  an  explicit  two-dimensional  finite-difference  technique. 
The  surface  response  procedure  solves  a  general  energy  balance  and  mass  bal¬ 
ance  based  on  a  film  coefficient  model.  The  ablation  events  are  computed 
through  consideration  of  the  complete  thermochemical  response  including  up  to 
four  specific  kinetically  controlled  reactions  with  a  carbonaceous  surface 
(all  other  possible  chemical  reactions  are  considered  to  occur  in  equilibrium) . 
This  approach  is  in  contrast  to  more  usual  correlation  schemes  or  "heat  of 
ablation"  approaches.  There  is  no  direct  feedback  coupling,  however,  between 
body  shape  change  and  boundary  condition  history?  that  is,  the  program  does 
not  include  routines  to  compute  pressure  distributions  and  convective  transfer 
coefficient  distributions  as  functions  of  body  shape  (these  are  input  a  priori 
as  functions  of  body  location  and  time) . 

The  basic  computer  program  is  called  ASTHMA,  for  Axi-Symmetric  Transient 
Heating  and  Material  Ablation."  The  specific  version  described  here  is  denoted 
as  Version  3  of  ASTHMA,  or  ASTHMA3. 

2.1.3  Coupling  to  Stress  Calculations 

The  program  does  not  include  thermal  stress  calculations  but  does  pro¬ 
vide  for  punched  card  output  which  may  be  used  directly  as  input  to  a  separate 
stress  program  (not  described  here) . 

2.1.4  Historical  Development 

The  ASTHMA 3  computer  code  described  here  constitutes  a  combination  of 
parts  of  other  programs.  Consequently  there  exists  a  large  amount  of  background 
literature  on  certain  parts  of  the  computer  code.  This  literature  will  be 
cited  below,  and  the  present  report  will  abbreviate  descriptions  of  aspects 
already  presented  elsewhere.  Consequently  the  reader  may  want  to  consult 
other  publications  for  more  detailed  expositions  of  those  aspects,  particularly 
the  surface  energy  balance  equations  (Refs.  6-10)  and  the  surface  thermochemical 
state  solution  procedure  (Refs.  1-5) . 
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2.2.1  Background  Remarks 

The  ASTHMA3  program  uses  a  fixed  finite  difference  heat  conduction  grid. 
Nodes  exposed  to  the  heated  surface  boundary  condition  may  erode  or  ablate. 

Nodal  boxes  at  this  surface  shrink  and,  as  they  approach  a  small  size,  are  dropped 
(amalgamated  into  the  next  node  down) . 

Schemes  of  this  type  must  be  formulated  with  some  care  to  avoid  oscillations 
in  predicted  surface  temperature.  Reference  11  describes  some  of  the  background 
of  the  development  of  the  node-dropping  scheme  used  in  ASTHMA,  which  has  been 
successful  in  avoiding  undesirable  oscillations. 

2.2.2  Nodal  Layout  and  Geometry 
2. 2. 2.1  General  Pattern 

The  geometric  shape  considered,  illustrated  in  Figure  2,  is  imagined  to 
be  divided  up  into  a  grid  pattern  in  the  customary  manner  for  finite  difference 
computation  schemes.  The  area  within  each  grid  "box"  will  be  termed  a  "node"; 
this  is  in  contrast  *-o  terminology  which  calls  each  corner  of  the  grid  network 
a  node.  The  thermal  capacity  of  each  node  is  imagined  to  be  lumped  at  a  single 
point  within  the  box;  this  point  will  be  termed  the  nodal,  center.  It  will  be 
convenient  to  assume  for  the  moment  that  the  nodal  center  may  be  located  any¬ 
where  within  the  nodal  box.  Actual  locations  of  the  nodal  centers  will  be  dis¬ 
cussed  later. 

For  convenience,  the  nodes  are  imagined  to  be  quadrilaterals,  so  that 
the  entire  nodal  network  is  an  assemblage  of  quadrilaterals.  For  bookkeeping 


Figure  2.  Sketch  of  Typical  Nodal  Layout 


2-3 


convenience,  the  network  is  imagined  to  be  "complete",  even  though  the  shape  of 
the  material  may  dictate  that  some  c.f  the  mesh  boxes  are  empty.  Both  nodes 
(boxes  or  centers)  and  mesh  corners  are  numbered  in  a  row  and  column  system 
which  will  be  labeled  as  an  m-n  system,  where  m  denotes  a  row  and  n  denotes  a 
column.  The  m-n  mesh  scheme  may  be  oriented  arbitrarily  with  respect  to  the 
physical  r-Z  coordinate  scheme,*  so  that  in  general  the  row-column  description 
might  seem  to  be  merely  a  descriptive  artifice;  however,  in  exploiting  the 
simplicities  associated  with  rocket  nozzle  geometries  it  has  been  assumed  in 
constructing  the  program  that  the  heated  surface  is  at  the  top  of  the  columns, 
as  indicated  in  Figure  2.  Thus  as  surface  recession  occurs,  the  boxes  at  the 
top  of  each  column  shrink;  the  other  boxes  remain  fixed.  This  limitation  that 
the  heated  surface  be  located  at  the  top  of  each  column  is  merely  a  convenience 
procedure  exploited  for  the  special  case  of  the  low  curvature  geometry  of  rocket 
nozzles.  For  high  curvature  bodies  this  restriction  becomes  inconvenient  and 
would  have  to  be  relaxed. 

The  sidewalls  and  the  backwall  are  either  insulated  or  in  communication 

through  a  simple  heat  transfer  coefficient  law  (plus  radiation)  to  a  "reservoir" 

at  T  .  Thus  the  boundary  condition  at  these  faces  do  not  involve  thermo- 
ires 

chemistry. 

Two  other  "convenience  limitations"  have  been  applied  to  the  layout  of 
the  nodal  grid.  First,  it  is  assumed  for  the  purposes  of  computing  thermal  con¬ 
ductance  between  nodal  centers  that  the  mesh  scheme  is  nearly  orthogonal,  so 
that  conductance  may  be  taken  as  conductivity  times  side-area  divided  by  length. 
Secondly  it  is  presumed  that  the  principle  direction  of  thermal  anisotropy  are 
aligned  with  the  nodal  mesh.  This  simplifies' computations  and  reduces  input 
requirements.  For  those  applications  in  which  the  heated  surface  intersects 
the  principal  directions  of  anisotropy  at  "difficult"  angles  this  restriction 
could  become  a  major  inconvenience  and  more  general  schemes  would  have  to  be 
divised  for  those  problems. 


it'—'"' 

That  is,  the  nodal  mesh  scheme  may  be  above  or  below  the  Z-axis  line,  may  be 
oriented  in  any  general  direction,  and  may  be  "bent"  or  shaped  in  any  manner 
convenient  to  the  user  (subject  to  the  limitations  described  in  the  paragraphs 
following  below) . 
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2 . 2 . 2 . 2  Geometry 

2. 2. 2. 2.1  Location  of  Nodal  Centers 

The  present  /ersion  of  ASTHMA  has  two  possible  nodal  location  schemes 
at  the  user's  option.  The  study  summarized  in  Reference  11  showed  that  to  mini¬ 
mize  surface  temperature  oscillations  during  node  dropping  operations  required 
that  no  thermal  capacity  be  associated  with  the  surface  temperature.  Conse¬ 
quently,  in  both  schemes  the  surface  temperature  has  no  thermal  capacity  and 
appears  as  an  extra  temperature  in  terms  of  the  number  of  nodes. 

In  the  back-shifted  nodal  location  scheme,  the  nodal  center  is  located 
at  the  center  of  the  "back  wall"  of  the  nodal  box,  that  is,  in  the  center  of 
that  side  parallel  to  the  heated  surface  and  farthest  from  it.  In  terms  of  the 
conventional  thermal  RC  network,  a  one-dimensional  (or  single  nodal  column)  rep¬ 
resentation  of  this  scheme  is  shown  in  Figure  3.  All  the  capacity  of  each  node 
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Figure  3.  Sketch  of  One-Dimensional  Thermal  RC  Network 

is  located  at  the  back  of  the  nodal  zone,  and  all  of  the  thermal  resistance  of 
this  zone  is  interposed  between  it  and  the  temperature  of  the  next  "higher" 

(i.e.,  m+1)  node.  Note  that  for  m  nodes  or  boxes  the  system  has  m+1  tempera¬ 
ture  points,  in  contract  to  most  schemes,  for  the  single  column  illustrated  in 
Figure  3. 

In  two  dimensions,  the  nodal  center  is  therefore  located  at  the  center 
of  the  "back  wall"  of  the  nodal  bjx,  that  is,  in  the  center  of  that  side  parallel 
to  the  heated  surface  and  farthest  from  it.  Thus  for  an  m  x  n  nodal  network 
schema,  there  will  be  m  x  n  nodal  temperatures  (associated  with  thermal  capac- 
ity)  plus  n  surface  temperatures  (not  associated  with  thermal  capacity) . 

Denoting  the  nodal  center  coordinates  as  rN  and  ZN,  the  coordinates 
for  a  back-shifted  node  m,  with  corner  coordinates  m,n  as  shown  in  Figure  4 
may  be  represented  in  terms  of  the  corner  coordinates  by  the  following: 
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_  rc,m,n  *  rc^min+l 


"N,m,n 


(1) 
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Figure  4.  Sketch  of  Nodal  Center  Location  for  Back-Shifted  Scheme 

The  second  scheme  employed  by  ASTHMA3  is  the  "centered"  scheme,  in  which 
nodes  are  located  in  the  arithmetric  center  of  the  nodal  box.  Figure  4  illus¬ 
trates  the  thermal  RC  network  for  a  one-dimensional  (or  single  nodal  column) 
version  of  this  scheme 
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Figure  5.  Sketch  of  One-Dimensional  Thermal  RC  Network; 
Centered  Scheme 


In  this  scheme,  the  nodal  point  of  the  surface  nodal  box  is  squeezed  toward  the 
back  wall  as  surface  recession  proceeds  so  as  always  to  remain  in  the  center  of 
whatever  nodal  volume  remains. 
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In  terms  of  the  nodal  box  corner  coordinates,  centered  nodes  are  given 


by: 

-  r°'lll'n  +  rc,m+l,n  +  rc,m+l,n+l  +  rc,m,n*fl 
rN,m,n  “  4 


Z  2  4*  Z  -}•  2 

„  _  c,m,n  c,m+l,n  c,m+l,n+l  c,m,n+l  , 

^N,m,n  “  4  w 

In  general,  the  centered  scheme  is  noticeably  more  accurate  than  the 
back*  shifted  scheme  (for  identical  node  sizes)  and  is  the  preferred  option. 

The  back-shifted  scheme  is  occasionally  useful  for  matching  nodal  point  to 
thermocouple  locations  or  to  interfaces  between  different  materials. 

2. 2. 2. 2. 2  Path  Lengths  for  Conductions 

In  this  and  the  following  section,  the  nodal  center  will  have  a  general. 

location  r„  zM  _  _  for  illustrative  purposes.  For  computing  thermal 
w,m,n  N,m,n 

conductances,  it  will  be  necessary  to  have  thermal  path  lengths  between  nodes. 
Since  material  properties  are  associated  with  each  nodal  box,  it  will  be  con¬ 
venient  first  to  consider  path  segments  inside  each  box.  In  general  there  are 
four  path  segments  of  .oterest,  as  shown  in  Figure  6.  The  program  computes  the 

lengths  L  „  and  L  _  in  the  m  direction  as  the  distances  between  the  nodal 
3  m,n,B  m , n , D 


Figure  6.  Illustration  of  Path  Lengths 
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center  and  the  centers  of  the  m+1  and  m  races  of  the  box.  (This  is  because  the 
nodal  center  is  always  on  the  line  joining  the  centers  of  these  two  faces.) 

Thus  the  paths  B  and  D  have  the  lengths 


m,n,B 


f  r  +  r 

c,m+l,n  c, m+1, n+1 

2  rN,m,n 


7i  +  2 

,  /  c,m+l,m  c ,m+l, n+1 
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r  , ,  +  r 
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2  rN,m,n 


it  +  z 

c,m,n+l  c,m,n  _  _ 

2  N,m,n 


-ii/ 2 


(6) 


The  paths  between  the  nodal  center  and  the  n  and  n+1  faces  (paths  A  and  C  in 
the  sketch)  are  assumed  to  end  at  points  on  those  faces  located  as  far  down  the 
face,  proportionally,  as  the  nodal  center  is  located  in  the  box.  In  the  back- 
shifted  scheme,  the  node  center  is  on  the  face  m,  for  which  the  paths  A  and  C 
end  at  the  corners  m,n  and  m,n+l,  respectively.  For  this  case,  we  have 


=  l(r  —  r  ) 2  +  (2,  -*Z  )^l 

,n,A  v  0,111,11  N,m,n  '  c,m,n  N, m,n  J 

,n,C  ”  [(rc,m, 


n+1  N,m,n 


)z  +  (Z 


—  z 

c,m,n+l  N,m,n 


-11/2 


(7) 


(8) 


In  the  arithmetical ly  centered  scheme,  the  paths  A  and  C  end  at  the  centers  of 
faces  m  and  m+1  respectively,  for  which  case  we  have 


m,n,A 


r  IT  +  IT 

c,m,n  c,m+l,n 
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z  +  3  x 
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Only  these  two  areas  need  to  be  computed  for  each  box,  since  the  areas  on  the 
other  sides  of  the  nodal  box  are  identical  with  the  areas  A  and  B  of  adjacent 
nodes. 


2. 2. 2. 2. 4  Volume 


The  volume  of  the  elemental  box  is  by  the  Second  Theorem  of  Pappus, 
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2. 2. 2. 2. 5  Geometric  Effects  of  Surface  Recession 


(13) 


For  nodes  adjacent  to  the  heat=d  surface,  side  B  may  move  due  to  surface 
recession  (ablation) ,  This  recession  reduces  the  thermal  resistance  between 
the  surface  point  and  the  adjacent  nodal  point,  increases  transverse  thermal 
resistances  (as  discussed  below) ,  and  reduces  the  thermal  capacity  associated 
with  the  nodal  center  of  the  nodal  box  adjacent  to  the  surface.  The  program 
assumes  that  the  surface  recession  occurs  so  as  to  maintain  the  ratios 


h  -  ^2  _  CJ2 

al  "  bl  C1 


(14) 


as  shown  in  Figure  8.  The  b  line  joins  the  centers  of  the  m+1  and  m  planes 
and  serves  to  define  the  locution  of  the  surface  point  S.  The  moving  corners 
m+l,n  and  m+1,  n+1  are  located  accordingly,  and  path  lengths,  areas  and  vol¬ 
umes  computed  as  before,  except  that  for  these  nodes  A  „  ,,  f  A  . 

m,n,A  m,n-l,c 
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Figure  8,  Sketch  of  Surface  Nodal  Box  Undergoing  Recession 
2 . 2 . 2 . 2 . 6  Surface  Shape 

As  noted  above  and  illustrated  in  Figure  8,  it  is  most  convenient  to 
consider  the  surface  points  on  the  heated  surface  as  being  always  located  on 
the  nodal  column  center  line,  that  is,  on  the  line  joining  the  center  points 
of  the  m+1  and  m  planes  ("parallel“  to  the  heated  surface) ,  where  m  is 
the  row  index  of  a  nodal  box  at  the  surface  and  m+1  is  the  local  corner  index 
of  the  heated  surface  (see  Fig.  8) .  The  location  of  the  surface  points  is 
all  the  information  needed  about  the  surface  Sor  those  ablation  problems  with 
a  specified,  input  surface  recession  rate  as  a  boundary  condition  (the  various- 
ablation  problem  boundary  condition  options  are  described  in  Section  2.3  below) 
since  for  those  problems,  it  is  most  convenient  to  specify,  as  input,  reces¬ 
sion  history  along  the  nodal  center  line.  The  nodal  grid  as  input  thus  serves 
to  define  the  various  nodal  column  center  lines,  and  the  recession  hl.tory 
for  each  column  then  defines  the  history  of  the  surface  points  in  a  perfectly 
straightforward  manner.  However,  ancthe;  important  heated  surface  boundary 
condition  option  involves  not  input  recessions  but  various  energy  and  chem¬ 
istry  information  sufficient  to  calculate  surface  mass  loss  from  energy  bal¬ 
ance  considerations  (Option  1,  discussed  xn  Section  2.3  below).  This  ablation 
calculation  does  not  produce  recession  rates  directly,  of  course?  instead 
it  produces  rates  of  mass  loss  from  the  surface.  It  will  prove  convenient, 
nevertheless,  to  adhere  to  the  concepv.  of  the  surface  point  which  moves  along 
the  column  center  line  as  recession  progresses.  To  define  the  surface  point 
motion  with  computed  mass  loss  rates  determined  from  this  general  energy- 
balance-determined  option  requires  knowledge  of  the  angle  between  the  local 
normal  to  the  surface  and  the  column  center  line,  since  computed  mass  loss 
may  be  translated  directly  into  recession  along  a  normal  to  the  surface.  Re¬ 
cession  along  the  normal  may  be  projected  into  the  nodal  column  center  line 
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once  the  angle  between  these  two  lines  is  known.  The  direction  of  the  sur¬ 
face  normal  may  conveniently  be  determined  from  the  slope  of  the  surface, 
that  is  the  surface  shape,  and  of  course  the  direction  of  the  column  center- 
line  is  known. 

It  is  obviously  not  safe  to  use  the  slope  of  the  heated  (top)  surface 
of  the  surface  nodal  box  to  obtain  the  surface  slope,  since  in  general  it  is 
neither  possible  nor  always  desirable  to  lay  out  a  nodal  grid  which  will 
"conform"  to  the  "real"  surface  for  the  entire  problem  history.  Therefore 
the  ASTHMA3  program  includes  both  linear  averaging  and  quadratic  curve-fit  sub¬ 
programs  which  compute  a  surface  shape  at  each  time  step  during  on*  solution 
after  examining  the  layout  of  the  surface  points. 

With  surface  slope  determined,  the  surface  movement  AS  computed  during 
the  time  step  may  be  projected  onto  the  nodal  box  centerline,  and  then  the 


Figure  9.  Sketch  of  Surface  Geometrical 
Relationships  (exaggerated) 

new  nodal  volume  computed,  as  indicated  in  the  exaggerated  sketch  of  Figure  9. 
Actual  surface  movement  during  a  time  step  is  limited  to  a  small  fraction  of 
the  nodal  thickness  to  ensure  a  good  approximation  to  the  conservation  of  mass. 

2. 2. 2. 2. 6.1  Surface  Shape  Curve  Using  Linear  Averaging  Technique 

If  the  linear  averaging  option  is  used,  referring  to  Figure  10,  the 


Figure  10.  Sketch  of  Surface  Points 
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slope  of  the  surface,  dr/dZ,  at  surface  point  n  as  the  average  of  the  slopes 
between -surface  points  n-1  and  n,  and  n  and  n-rl.  Thus 


r\  _  1  /rn  ~  rn-l  ,  rn*l  ~  rn\ 
Z)n  "  2  \Zn  "  Vl  '  Zn+1  "  Zn) 


Note  that  n-1,  n,  and  n+1  are  surface  points,  not  nodal  corner  points.  (Prob¬ 
lems  with  only  one  nodal  column  have  only  one  surface  point.  This  requires  the 
program  to  abandon  the  average  slope  scheme  and  to  use  the  nodal  box  heated  sur¬ 
face  slope  as  the  surface  slope.) 

2. 2. 2. 2. 6. 2  Surface  Shape  Curve  Using  Quadratic  Technique 

The  quadratic  curve  fit  slope  averaging  routine  known  as  SL0PQ  is  de¬ 
scribed  in  Reference  13.  This  routine  gives  better  results  for  relatively 
smooth  and  flat  (small  dr/dZ)  shapes  common  to  rocket  nozzles;  it  usually  per¬ 
forms  poorly  if  surface  points  are  closely  spaced  in  the  Z  direction  (large 
dr/dZ)  and  gives  meaningless  answers  for  double  valued  r(Z)  shapes.  For  these 
latter  shapes,  the  linear  fit  routine  should  be  used. 

2.2.3  Internal  Conduction  Parameters 

The  thermal  resistances  between  nodes  and  the  thermal  capacity  of  the 
nodes  are  calculated  each  time  interval  from  the  material  properties  of  the 
node  corresponding  to  its  temperature  at  that  time.  The  material  properties 
density  (p) ,  specific  heat  (c) ,  conductivity  (k) ,  and  emissivity  (e)  are  in¬ 
put  as  table  look  up  functions  of  temperature.  Linear  interpolation  is  em¬ 
ployed  for  material  property  determination  at  temperatures  intermediate  to 
those  tabulated.  Constant  thermal  contact  resistances  may  be  specified  be¬ 
tween  any  or  all  nodes.  The  nodal  capacities  and  resistances  are  calculated 
as  follows: 


C  .  =  p  .c  „V 

m,n,0  ^m,n,6  m,n,0  m,n 


^  i _  — m f_n f C  +  _  m, n+l,A  +  R* 

m,n, A,  0  A»n,n+1,A  |km,n,0  km,n+l,0  m,n,B 


+  Lm+1, n, D  +  R* 


Rm,n,B,0  Am+i,n,B  \km,n,0  km+l,n,0  +  Rm'n'A 
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The 


Figure  12  shows  the  locations  of  resistances  R  _  .  and  R  _  Q 
3  ra,n,A»0  m,n,B,0 


m,n-l 


•  m-1 , n 


Figure  12.  Sketch  of  Thermal  Resistance  Nomenclature 

resistance  on  the  other  two  sides  of  the  nodal  box  are  calculated  when  the 
quantities  for  the  adjacent  nodes  are  calculated.  For  nodes  adjacent  to 
the  surface,  however,  A^  n  c  r  A^  n+1  A  generally  (refer  to  Fig.  7  for  area 
nomenclature) ;  for  these  nodes 

R  _  _Lm,n,_C _  +  Lm,n+1,  C  +  Rm,nxB  (19) 

m,n,A  Am,n,Ckm,n,0  Am,n+1, Akm,n+1, 6  Am,n, C 

It  should  be  noted  that  anisotropic  thermal  conductivity  values  k  and  k'  are 
used  in  Equations  (17)  and  (18) ,  respectively. 

A  subroutine  has  been  added  to  ASTHMA3  which  decreases  inaccuracies  in 
the  calculated  conductances  due  to  non-orthogonality  problems.  Figure  13 
shows  the  general  situation  in  a  conductance  calculation  (selecting  the  "along 
rows"  direction  as  an  example) .  » 


Figure  13.  Sketch  of  Conduction  Quantities 
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The  conductance  used  by  the  program  is 


U  = 


A 


km,n  km,n+l 


An  obvious  improvement  in  the  number  is  a  modified  conductance 


(20) 


U'  =  U  sin  6 


(21) 


This  modified  conductance  accounts  exactly  for  some  limited  non-orthogonality 
problems:  those  in  which  the  basic  nodal  point  array  is  orthogonal  but,  for  one 
reason  or  another,  the  nodal  box  sides  are  not  orthogonal  to  the  point  array. 
For  funadmental  non-orthogonalities  of  the  nodal  point  array,  the  conductance 
improvement  does  increase  the  accuracy  of  computed  heat  fluxes  by  a  factor  of 
roughly  1/sin  0.  Fundamentally,  however,  this  situation  requires  a  nine  point 
difference  scheme  not  treated  by  ASTHMA. 

2.2.4  In-Depth  Conduction  Solution 

The  in-depth  conduction  solution  is  the  explicit  finite  difference  type 

often  employed  for  transient  heat  conduction  analysis.  The  temperature  of 

node  m,n  at  time  0 1 (T  .,)  is  obtained  by  application  of  the  finite  differ- 

iri/  n  *  y 

ence  energy  balance  and  rate  equations  to  the  nodal  volume. 

Solving  for  T  « , ,  one  obtains: 


T 


m,n, 0 1 


T  T  T  T 

ttH-l,n,  8  +  m,n+I, 0  +  m-l,n,  9 _  +  mtn-l,  8  _  T 

Rm,n,B,0  Rm,n,A,0  Rm-l,n,B,8  Rm,n-1,A,0  m,n,e 


Rm,n-1,A,8  Rm,n,B, 0  Rm,n,A,0  Rm-l,n,B,0 


c77:  +  Tm,„,8  <22> 

m*  n  /  0 


In  the  program  this  equation  is  used  to  obtain  "new"  temperatures  for 
all  nodes  except  those  adjacent  to  the  heated  surface  and  for  back  wall  nodes. 
Nodes  adjacent  to  the  heated  surface  are  linked  to  the  surface  temperature 
implicitly  and  hence  a  special  procedure  is  used  for  the  temperature  of  these 
nodes,  as  described  in  the  next  section. 

Back  wall  nodes  include  a  quantity 


k^bw^res  Tm,n,0^  +  0ebwAbw ^Tm,n, 0  "  Tres^ 
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inside  the  brackets  of  Equation  (22)  .  Any  node  other  than  a  heated  surface 
node  may  be  called  out  as  a  back  wall  node,  and  A^w  may,  with  appropriate  input 
choices,  be  assigned  as  one  of  the  two  side  areas  Am  „  and  A_,  or  as  the 
back  area  A^  n  D«  Physically,  such  assignments  only  make  sense  if  the  back  wall 
node  thus  specified  is  located  at  the  sides  or  back  wall  of  the  tot«.l  material 
considered,  or  adjacent  to  a  void  space  within  the  material. 

The  explicit  relation  (22)  imposes  the  familiar  stability  restriction  on 
the  time  step  size  AG  =  6*  -  6.  The  ASTHMA3  prograi  automatically  employs  a  con¬ 
servative  stability  equation  for  interior  nodes: 


AG 


m,n 


R 


R 


m,n,A  m,n,B 


R_ 


m,n-l,A  m-l,n,B 


(23) 


Normally  for  stability,  the  input  parameter  n  is  less  than  unity.  Surface  nodes 

are  not  considered  for  time  interval  calculation,  as  will  be  explained  below; 

back  wall  nodes  include  the  terms  A.  (h/2  +  4oe,  T3  Q)  in  the  denominator. 

ow  bw  m ,  n ,  G 

The  automatic  stability  criterion  calculation  may  be  suppressed  for  any 
node  if  the  user  is  sure  that  the  allowed  time  step  for  that  node  will  never  be 
the  minimum  one  for  the  system.  This  saves  some  computer  time. 

Alternatively,  the  stability  criterion  calculation  can  be  suppressed 
entirely.  If  it  is  not  used,  then  AG  must  be  specified.  It  is  used,  then  any 
or  all  nodes  in  the  conduction  network  may  be  specified  for  time  interval  de¬ 
termination. 

2.2.5  Temperatures  of  Surface  Nodes  and  Surface  Points 

Temperatures  of  surface  points  are  determined  either  by  assignment  (Op¬ 
tion  2)  or  by  the  surface  energy  balance  described  in  Section  2.3  (Options  1 
and  3) .  The  surface  energy  balance  determines  the  new  surface  temperature  of 

the  nth  column  Tl  with  an  implicit  iteration  technique.  Stability  considera- 
Bn 

tions  dictate  that  the  first  node  temperature  also  be  treated  implicitly,  and 
that  any  transverse  heat  conduction  link  (across  columns)  for  surface  tempera¬ 
tures  must  be  implicit.  This  latter  requirement  is  a  complex  one  to  meet  as 
columns  recede;  hence  the  surface  temperature  points  are  not  linked  transversely. 
Figure  14  shows  the  implicit  and  explicit  heat  conduction  paths  for  two  typical 
columns. 
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Figure  14.  Sketch  of  Implicit  and  Explicit  Temperature  Links  in  Finite 
Difference  Solution,  for  Two  Typical  Situations 

For  this  scheme,  the  equation  for  the  new  temperature  of  a  surface  node 
(not  a  surface  point)  is 


T 


m,n,9 ' 


Tm,n+1,9  -  J'm,n,8  Tm,n~l,8  ~  Tm,n,9 

R  » 

m,n,A, 0  m,  n-1,  A,  9 


Tm-l,n,8  ~  Tm,n,8  +  ^s.n^'  ~  ^nun.g* 


A0 


_ - 1-  T_  a 

C  *  m, n, 0 
m,  n ,  0 

(24) 


This  is  linked  to  the  surface  energy  balance  through  the  surface  point  temper¬ 
ature  Ta  „ , .  Formally  we  have  a  relation  between  the  two  unknowns  as 
o ,  n ,  0 
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T  '  =  f  (T  . ) 

m,n,0  v  S,n,0‘' 


(25) 


The  surface  energy  balance  has  the  general  form  (aj  described  in  the  next 
section) 


convection  and  chemical  energy  terms  (m  . , ) 

s ,  n,  u 

+  radiation  to  wall  -  radiation  out  (Tc  , ) 

o  f  n  f  u 


1S,n, 0  1  Vn.91 

Rm,n,B,0 


f  2 ^m,n,  0  1  TS,n, 8 


(26) 


where  the  parentheses  denote  functional  relationship.  Relation  (25)  may  be 

substituted  into  Equation  (26)  to  give  the  non-linear  surface  energy  balance 

Equation  for  Tc  Q , .  When  Tc  a ,  has  been  found  from  this  equation,  the 

new  surface  node  temperature  Tm  n  g ,  may  be  found  from  Equation  (25).  For 

surface  energy  balance  options,  the  method  for  finding  T„  is  described 

o  /  ri|  o 

in  Section  2.3  below.  In  Option  2,  Tg  n  Q,  is  known  immediately  and  Equation 

(24)  then  determines  T  „  a ,  at  once. 

m,n,  0 

In  some  rare  cases  with  very  high  recession  rates,  the  implicit  linkage 
described  above  is  not  sufficiently  tight  to  prevent  destructive  oscillations. 
Therefore  this  linkage  was  supplemented  with  a  second  linkage  which  implicitizes 
the  first  node  temperature  in  the  conduction  path  from  the  first  node  to  the 
next  deeper  node. 

In  Equation  (24),  this  second  option  has  0'  in  place  of  0  on  Tm  n  0  in 

the  n  b  e  term*  Similarly,  the  nodal  equation  for  node  m-l,n,0  has  0'  in 

place  of  0  on  T  «  (i.e.,  this  heat  conduction  rate  must  be  identical  in  the 

two  nodal  energy  balances).  In  this  scheme,  the  determination  of  T 1  .  is 

m — /  n  /  v 

delayed  until  Te  „  a '  and  T  „  Dl  are  determined  as  described  above. 

J  S,n,o  m,n , e ' 

This  linkage  option  has  always  performed  stably.  It  is  slightly  slower 
than  the  first  option. 

2.3  SURFACE  BOUNDARY  CONDITION  ASPECTS 

The  ASTHMA  program  ca:  account  for  three  different  types  of  heated  sur¬ 
face  boundary  conditions: 

1.  Surface  energy  balance  determined  recession  rate,  employing  general 
film  coefficient  model  for  boundary  layer  transport  of  energy  and 
mass,  and  very  general  thermochemical  relations  at  the  surface 
(Option  1) . 
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2.  Specif if ed  surface  temperature  and  recession  rate  histories  (with  a 
different  history  at  each  surface  point  of  interest)  (Option  2) . 

3.  Simplified  energy  balance  to  determine  surface  temperature  including 
only  radiation  terms,  no  recession  allowed  (Option  3  -  "cooldown 
option") . 

Option  2  is  so  simple  that  it  requires  no  discussion  other  than  that 
given  in  the  user's  manual  below.  Option  1  is  complex  and  has  many  interest¬ 
ing  aspects,  some  of  which  are  discussed  below.  A  brief  discussion  of  Option  3 
follows. 

2.3.1  General  Aspects 

The  convective  surface  energy  balance  option  (Option  1)  of  the  program 
performs  the  energy  balance  illustrated  in  the  sketch  of  Figure  15.  The  energy 
balance  is  performed  for  the  indicated  control  volume  fixed  to  the  receding 


qabs/AB 


Figure  15.  Representation  of  Surface  Energy  Terms  During  Ablation 

surface.  Energy  fluxes  leaving  the  control  volume  include  conduction  into  the 
material,  radiation  away  from  the  surface,  energy  in  any  flow  of  condensed  phase 
material  such  as  mechanical  removal,*  and  gross  blowing  at  the  surface.  Energy 
inputs  to  the  control  volume  include  radiation  in  from  the  boundary  layer  and 
enthalpy  fc-.ix  due  to  the  convection  of  material  to  the  surface  associated  with 
surface  recession.  The  special  flux  qa£S  represents  an  energy  addition  not 
proportional  to  surface  area,  such  as  a  narrow  particle  or  radiation  beam  or 
induction  heating.  The  final  input  in  u.he  sketch  is  denoted  q^ff*  It  includes 
all  diffusive  energy  fluxes  from  the  gas  phase  boundary  layer  and  the  form  of  a 
correlation  equation  associated  with  the  convective  film  coefficient  model. 


The  energy  balance  procedure  has  indeed  been  constructed  to  include  mechanical 
terms  but  in  this  version  of  ASTHMA  failing  is  not  allowed  (in  order  to  exploit 
input  conveniences)  since  it  is  not  provided  for  in  the  surface  thermochemistry 
code  GASKET  associated  with  this  version  of  ASTHMA. 
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The  computation  of  the  surface  energy  balance  requires  from  the  in-depth 
solution  a  relation  between  the  surface  temperature  and  the  rate  of  energy 
conduction  into  the  material,  <3ccn(j*  This  relation  derives  naturally  from 
the  finite  difference  energy  balance  for  the  node  just  under  the  surface.  With 
this  information  the  surface  energy  balance  considerations  allow  determination 
of  the  thermochemical  erosion  rate  m  and  surface  temperature  T^.  It  will 
be  useful  to  keep  in  mind  that,  from  this  point  of  view,  the  purpose  of  the 
in-depth  solution  at  any  instant  is  to  provide  information  about  qcon(3  (Tw)  . 

The  surface  energy  balance  equation  may  be  written  as  (ignoring  the  term  q  ^ , 
which  is  only  rarely  of  interest) 


qdiff  +  qrad  +  ™chc  -  "  (Pv)whw  “  £  %hi  "  qcond  =  0 


where 


(pv)  -  rn  - 
y  w  c 


Y.  v 


The  relation  qconcj  =  f(Tw)  is  delivered  by  the  in-depth  solution, 
pendencies  of  interest  are 


Other  de- 


hc  -  hc(Tw> 


=  q  .  (t  j 

^rad'  w 
out 


For  the  other  terms,  we  may  write  in  general 


V  qdiff  '  qrad  '  ^w  '  ^  ,  mr 

in  i  * 


functions  of  boundary-layer-  (31) 
edge  enthalpy,  pressure,  up¬ 
stream  events,  laws  for  con¬ 
servation  of  chemical  elements, 
chemical  equilibria, and/or 
kinetic  relations,  m  ,  local 
boundary  layer  aerodynamics 


Relations  of  the  type  of  Equation  (31)  come  in  many  forms,  and  may  even 
take  the  form  of  exact  solutions  to  boundary  layer  and  surface  chemical  state 
routines.  In  the  present  program,  the  relations  (31)  are  computed  by  the  sepa¬ 
rate  chemistry  program  (ACE  or  GASKET) ,  for  use  by  Equation  (27)  in  ASTHMA3 
based  on  a  film  coefficient  model  of  the  boundary  layer  transport  events.  The 
ACE  program  is  fully  described  in  Reference  4,  and  GASKET  is  described  in 
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Reference  5.  Since  relations  (31)  are  obtained  by  film  coefficient  based  pro¬ 
grams  (ACE  and  GASKET) ,  the  main  energy  balance  (27)  used  by  the  ASTHMA  program 
is  also  of  a  film  coefficient  type. 


When  the  boundary  layer  transport  aspects  of  the  problem  are  modeled  by 
a  film  coefficient  scheme,  then  both  Equations  (27)  and  (31)  can  be  normalized 
on  the  mass  transfer  coefficient  in  the  customary  manner.  Equation  (27)  becomes 


ydiff 


qrad  **  qrad 

in _ out 

P  u  C., 

Ke  e  M 


B'h 
c  c 


-  B'h 


w 


rn 


■peueCM 


qcond 

peueCM 


=■•  0 


(32) 


And  Equation  (31)  takes  the  general  form 


qrad  r — a  m 

qdiff  in  h  \  l  h 
w'  p  u  C  '  p  u  c  '  w  '  /  ,  p  u  C„ 


e  e  M  e  M 


e  e  M 


functions  of  bound¬ 
ary  layer  edge  en¬ 
thalpy,  pressure, 
laws  for  conservation 
of  chemical  elements, 
chemical  equilibrium 
and/or  kinetic  rela¬ 
tions. 


(33) 


The  generation  of  the  Equation  (33)  relationships  is  the  goal  of  the 
surface  thermochemistry  code  (ACE  or  GASKET)  as  noted  above.  For  the  present 
it  may  be  observed  that  the  energy  balance  coupling  to  the  in-depth  solution 
(ASTHMA3)  for  each  time  step  proceeds  as  follows:  an  initial  guess  of  the  dimen¬ 
sionless  mass  removal  rate  B1  is  obtained  in  some  manner.  With  this  B',  the 

c  .  c 

quantities  <5(jiff/peueCM'  hw'  and  Tw  are  obtained  from  the  sur“ 

face  thermochemistry  solution.  The  quantities  hc  and  qrad  out  are  then  formu¬ 
lated  using  the  Tw  so  obtained.  The  surface  energy  balance  is  then  computed, 
the  qcQnd  as  a  function  of  Tw  having  been  provided  by  the  in-depth  solution. 

In  general,  however,  the  sum  of  the  terms  will  not  equal  zero  but  some  error. 

An  iteration  procedure  is  then  used  to  select  successively  better  estimates 
of  B'  which  drive  the  error  to  zero.  Experience  shows  that  Newton's  procedure, 
in  which  the  derivative  of  the  error  with  respect  to  is  used  to  compute  the 
next  guess  for  B',  gives  good  results. 


2.3.2  Computational  Approach  to  Convective  Energy  Balance 
2. 3. 2.1  General  Description  of  Approach 

It  is  evident  that  each  iteration  in  the  search  for  a  surface  energy 
balance,  if  performed  as  described  above,  would  require  a  new  surface  chemistry 
solution,  generally  in  the  near  neighborhood  of  many  such  previous  solutions. 

This  suggests  that  a  tabular  approach  in  which  the  surface  state  solutions  are 
done  only  for  preassigned  values  would  offer  significant  computational  econo¬ 
mies,  since  values  of  dependent  quantities  at  values  between  table  points 
would  be  obtained  by  interpolation  rather  than  by  additional  thermochemistry 
solutions.  Referring  to  Equation  (33),  we  see  that  such  a  table  would  have 
three  independent  variables.  Two  obvious  ones  are  pressure  and  the  ablation 
rate  B^.  The  non-obvious  third  independent  variable  involves  chemical  kinetic 
control.  If  surface  reactions  are  kinetically  controlled,  then  the  discussions 
of  the  ACE  and  GASKET  user's  manuals  (References  4  and  5  respectively)  demon¬ 
strate  that  the  kinetically  controlled  relations  show  up  in  the  array  of  chemis¬ 
try  relations  to  be  solved  in  a  very  particular  form:  the  pre-exponential  fac¬ 
tors  on  forward  rate  coefficients  appear  divided  by  the  transfer  coefficient 
peueCM*  Therefore  the  thermochemistry  solution  routine  must  be  provided  either 
with  both  numners  (pre-exponential  factors  and  transfer  coefficient)  or  alterna¬ 
tively  with  their  ratio. 

The  various  Aerotherm  thermochemistry  cedes  differ  at  this  point,  and 
succeeding  discussions  will  be  specialized  to  the  GASKET  code,  which  has  rate 
coefficients  built  into  the  data  of  the  code  and  requires  the  transfer  coeffi¬ 
cient  PeueCM  to  be  input.  Thus  the  only  kinetics  parameter  evident  to  the  user 
is  the  transfer  coefficient  PeuecM/  and  this  will  be  the  third  parameter  in  the 
table. 

The  course  of  a  typical  ASTHMA  solution  history  at  one  axial  station  might 
be  represented  as  shown  in  Figure  16.  As  time  proceeds,  solutions  progress 
through  the  space  of  the  table  independent  variables  B^ ,  PeueCM'  and  The  dots 
in  the  picture  represent  solutions  satisfying  the  surface  energy  balance  as  time 
progresses.  Surrounding  these  points  are  a  number  of  other  points  examined  dur¬ 
ing  the  iterations  to  satisfy  the  surface  energy  balance.  For  any  iteration,  the 
solution  procedure  finds  itself,  so  to  speak,  within  a  cube  formed  by  the  bracket¬ 
ing  tabular  values  of  ,  Peuecj^/  and  The  dependent  quantities  Tw, 
peueCM'  an<^  £mrV12/peueCM  ^ave  been  precalculated  for  these  tabular  points; 
relevant  values  of  these  quantities  for  the  current  iteration  values  of  B^, 
peueCM'  an<^  P  can  then  be  formed  by  interpolation  inside  the  cube  and  the  sur¬ 
face  energy  balance  equation  calculated.  If  the  energy  balance  is  not  satisfied 
to  some  preselected  degree  of  accuracy,  a  new  value  of  B^  can  be  selected  and 
the  process  repeated. 
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Figure  16.  Representation  of  Course  of  Independent 
Variables  in  Surface  History  at  One 
Axial  Station. 


The  tabular  approach  was  selected  in  general  because  the  interpolation 
feature  drastically  reduces  the  number  of  surface  state  calculations  while  at 
the  same  time  allowing  sufficient  accuracy.*  The  precalculated  table  approach** 
was  chosen  for  the  following  reasons: 

1.  In  parametric  studies,  tables  once  generated  are  usable  for  many  dif¬ 
ferent  problems,  yielding  even  greater  economy. 

2.  For  many  problems,  it  is  difficult  to  specify  a  priori  an  "adequate" 
array  of  independent  tabular  values  B^,  peuoCj,  and  P.  An  examina¬ 
tion  of  the  precalculated  surface  tables  before  execution  of  the 
actual  ablation-problem-plus-in-depth-solution  can  reveal  if  there 
are  any  "holes"  in  the  tables  in  areas  where  energy  terms  are  vary¬ 
ing  rapidly.  Desirable  table  points  can  be  added  before  the  .n-depth 
response  run. 

3.  The  surface  tables  are  frequently  of  independent  interest  ii  them¬ 
selves  for  judging  the  ablation  effects  under  various  conditions.. 

4.  Finally,  without  the  precalculated  tabular  approach,  the  occasional 
nonconvergent  surface  chemistry  solution  would  stop  the  entire  in- 
depth  solution  process.  With  the  precalculated  table  approach,  such 
solutions  are  automatically  weeded  out  of  the  tables  without  damage 
to  the  subsequent  in-depth  solution. 

On  the  other  side  of  the  ledger,  some  disadvantages  in  the  precalculated 
table  approach  are  evident  however: 

1.  Figure  16,  which  is  a  realistic  schematic  view  of  a  typical  calcula¬ 
tion  history,  indicates  that  many  of  the  laboriously  calculated  and 
assembled  surface  state  points  in  the  table  are  never  used  in  the 
course  of  a  given  solution. 

2.  The  "mechanical"  linkage  between  the  surface  state  solution  and  the 
in-depth  solution,  i.e.,  the  transfer  of  punched  card  surface  state 
output  to  input  of  the  in-depth  program,  leads  to  computing  delays 


For  example,  a  typical  2000  time  step  problem  with  the  usual  5  iterations  per 
time  step  would  require  10,000  surface  state  calculations,  which  require  about 
1  second  each  for  machines  in  the  1108  speed  class.  A  single  pressure  table, 
on  the  other  hand,  involves  only  about  200  state  solutions  (20  char  rates 
times  10  transfer  coefficient  values),  a  factor  of  50  improvement.  Multiple 
pressure  tables  reduce  this  advantage,  but  usually  only  a  few  pressures  are 
required. 

**As  opposed  to  a  direct  coupling  scheme  in  which  tabular  dependent  entries  were 
only  calculated  as  needed. 
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and  occasionally  to  cross  input  blunders  (wrong  decks,  missing  parts 
of  decks,  etc.). 

3.  Significant  storage  is  required  by  the  large  precalculated  table. 

These  disadvantages  are  usually  outweighed  by  the  advantages  of  the  pre¬ 
calculated  table  approach,  however.  Hence  this  same  well  tested  approach  has 
been  adopted  for  the  ASTHMA3  program. 

2. 3. 2. 2  Pressure  Interpolation  Not  Needed 

The  preceding  general  discussion  has  described  how  interpolation  in  a 
precalculated  table  with  three  independent  variables  can  be  carried  out.  Since 
the  ASTHMA3  code  described  here  was  designed  primarily  for  rocket  nozzle  applica¬ 
tion,  a  useful  simplification  is  possible.  In  rocket  nozzles,  pressure  rarely 
varies  by  an  appreciable  amount  with  time  at  a  given  axial  station.  Therefore 
surface  energy  balance  solutions  need  not  interpolate  on  pressure,  even  though 
it  is  necessary  to  have  on  hand  several  thermochemistry  solutions  at  different 
pressures.  It  will  be  sufficient  for  each  axial  station  to  refer  to  the  appro¬ 
priate  pressure. 

This  simplification  allows  more  economical  interpolations  {two-dimensional 
instead  of  three-dimensional) .  It  is  exploited  in  ASTHMA  by  the  assignment  of 
each  surface  node  to  a  given  pressure  table  by  simple  numerical  flag  which  re¬ 
fers  to  the  number  of  pressure  tables  counced  in  the  order  in  which  they  are 
read  in  (see  Section  3.1.6  below). 

2.3.3  Forms  of  the  Film  Coefficient  Model  Surface  Energy  Balance 
Equation  Used  in  Option  1 

In  making  the  surface  energy  balance  calculation  described  in  Section 

2.3.2  above,  the  ASTHMA3  program  uses  a  variety  of  surface  energy  balance  equa¬ 
tions  of  the  general  form  of  Equation  (32).  The  basic  equation,  derived  in 
Reference  (6)  with  analogy  arguments,  is* 

p  u  Cr..(H  -  h  )  .  +  p  u  <1, 

Ke  e  h  r  w  edge  1  e  eH 

gas 

+  ^abs/AB  +  V*rad  "  FoeTw  "  ^cond  =  0  <34> 

It  should  be  noted  that,  like  all  film  coefficient  expressions,  this  equation 
is  not  a  universally  valid  one.  It  is  well  established  for  frozen  boundary 


A  • 

ASTHMA 3  assumes  m  =0.  This  excludes  use  of  the  "fail  temperature"  option  in 
ACE  and  GASKET  whenZ  preparing  ASTHMA3  input. 


V(Ki  -  *.  ) 

1—lQ  W 


h. 


w 


B'h  -  B'h 
c  c  c  w 


layer,  catalytic  wall  problems  with  no  net  mass  transfer  (Refs.  13,  14)  . 

It  is  presumed  accurate  for  reactive  boundary  layers  as  well  (see  a  discussion 
of  similar  considerations  in  Ref.  15)  -  For  problems  with  net  mass  transfer 
(ablation)  ,  which  are  of  interest  here,  the  equation  is  less  well  established, 
but 

1-  Has  been  derived  with  “respectable”  analogy  arguments  (Ref.  6) 

2.  Compares  well  with  similar  results  derived  by  Lees  (Ref.  14) 

3.  Has  given  good  predictions  when  compared  to  some  experimental  data 
and  to  results  of  "exact"  boundary  layer  solutions. 

4.  Gives  results  independent  of  the  choice  of  enthalpy  datum  state 
(an  important  criterion) 

5.  Reduces,  for  C,.  =  c„,  to  a  widely  accented  energy  balance  law 

ri  n  ~ 

(Refs.  16,  17)  for  simultaneous  heat  and  mass  transfer 


Equation  (34)  has  been  generalized  in  Reference  7  to  the  case  of  unequal  mass 
diffusion  coefficients,  again  by  analogy  arguments.  The  result  is 


p  uC_(H  -  h  )  .  +  p  -u  C„ 

e  e  H  r  w'edge  we  e  M 

gas 


E(Z*  -  Z*  )h^w  +  B'h  -  B*h 

1e  xw  1  c  c  c  w 


l  x 


where 


+  <WAB  +  “w^rad  -  FoeTw  ^cond  =  °  05) 
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(36) 


and 
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(37) 


(38) 
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The  factors  F-  in  these  eqaafcioos  derive  from  the  particular  relation  between 
the  binary  diffusion  coefficients  which  nnst  hold  if  the  governing  differen¬ 
tial  equations  are  to  reduce  to  the  fares  frcta  which  Relation  (32*  can  be 
inferred.  Shis  relation  is 


•^ij 


F-F- 
-  3 


(39) 


and  can  be  regarded  as  an  accurate  correlation  of  experimental  data  for  the 
binary  diffusion  coefficient  5—.  She  quantity  D  is  a  constant  for  a 
given  pressure.  She  constants  F-  depend  weakly  on  temperature. 

Equation  (35)  is  somewhat  less  well  founded  than  the  corresponding  Equa¬ 
tion  (34)  but  does  have  plausibility.  Furthermore,  the  equation  reduces  to 
Equation  (34)  for  equal  diffusion  coefficients,  as  it  should,  and  Equation 
(35)  can  be  shown  to  be  independent  of  the  enthalpy  datum  state  and  thus 
fulfills  a  basic  physical  requirement. 

The  GASKET  program  provides  to  the  ASTHMA 3  program  the  four  dependent 
quantities  (Tw,P) ,  ^(B^P) ,  %  Z±eh*w  (B;,P)  and  £  zf^h^O^P) 


S)  stored  in  one 


for  a  given  value  of  kj..  With  the  quantity  h~ . 

!■  cGyc  OdS  W* 

table  with  independent  variables  Tw  and  p,  and  the  other  three  dependent  quanti¬ 
ties  stored  in  a  table  with  independent  variables  9’,  put,  and  P,  the  ASTHMA3 
program  may  find  surface  energy  balances  at  any  surface  point  and  time  (i.e. , 
given  pressure  and  PeuecM)  by  (1)  selecting  the  correct  P  table,  (2)  noting 
Peue0’H,  (3)  finding  the  dependent  quantities  in  the  two  tables  at  peueCj,  values 
bracketing  the  actual  PeueCH  and  for  same  B^,  (4)  interpolating  on  Peuec|,j»  (5) 
forming  the  energy  balance  equation  (35)  and  noting  the  error  or  departure  from 

The  Newton-Raphson  method  is  used 

for  selecting  the  next  guess  of  B'  in  step  (6) .  When  that  B’  yielding  an  energy 

c  c 

balance  is  determined,  the  new  Tw  is  known  as  a  dependent  quantity;  thus  the  new 
surface  temperature  is  determined. 


zero,  (6)  selecting  a  new  B^  for  another  try 


2.3.4  .input  and  Correction  of  Heat  Transfer  Coefficient 

To  employ  the  film  coefficient  formulation  just  described,  the  program 
user  must  provide  the  program  with  values  of  the  heat  transfer  coefficient 
°eueCH  as  functi°ns  °f  time.  Two  practical  problems  must,  be  settled  in  this 
respect: 

1.  How  is  C 


M 


related  to  C..? 


2.  Can  both  and  be  specified  as  functions  of  edge  conditions 

(i.e.,  of  time)  independent  of  the  subsequent  problem  solution 
(i.e.,  mass  transfer  rates  and  body  shape)? 
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Tn>  answer  to  the  first  question  it  nay  be  stated  that  within  the  present 
fcrmmiaticn  it  is  adequate  to  take  the  ratio  C^/fL.  as  a  constant.  She  value 
of  thia.  constant  is  a  measure  of  the  ratio  of  the  mean  mass  transfer  aspects 
of  the  boundary  layer  to  the  mean  heat  transfer  aspects.  For  equal  mass  dif¬ 
fuse, -»n  coefficients,  a  vast  amount  of  experimental  data  (as  summarized,  for 
example,  in  Ref.  17  suggest  the  correlation  C^/C..  =  Le^.  it  cay  be  hvpofche- 
sjLved  that  far  unequal  mass  diffusion  coefficients  the  same  procedure  may  be 
employed  v-ith  the  Lewis  Eumber,  Le,  defined  by  the  procedure  set  forth  in 
Reference  7  involving  D.  Stas,  the  input  to  the  program  consists  of  a  tine 
table  of  values  for  p^u^C..  and  the  constant  factor  C^./C„- 

Za  answer  to  question  2.,  changes  of  CL.  with  body  shape  are  not  usually 
of  interest  for  rocket  nozzle  problems  and  hence  are  not  accounted  for  in  the 
A33SHA  program.  A  core  important  problem  concerns  the  dependence  of  C__  on 
the  actual  rate  of  mass  transfer.  This  problem  has  been  ignored  up  to  now, 
the  implication  being  that  C„  is  determined  bv  the  bound arv  layer  edge 
aerodynamics  alone.  This  is  known  to  be  incorrect.  The  value  of  C„  depends 

*  O 

fairly  strongly  on  n.  If  we  denote  the  with  m  =  0  as  cHq»  this  depen¬ 

dence  is  shown  by  both  data  and  analysis  to  be  accurately  represented  by 


where 


cp  ^ 


2Xm 


peueCH, 


(40) 


(41) 


This  correction  is  built  into  the  program. 

Most  user's  prefer  a  value  of  0.4  for  X  _n  turbulent  flow,  and  a  value 
of  0.5  for  X  in  laminar  flow.  These  values  appear  to  be  acceptably  accurate  for 
B'  values  below  about  2  provided  the  molecular  weight  of  the  gas  at  the  wall  is 
at  least  roughly  equal  to  the  molecular  weight  of  the  environment  gas.  Very  high 
B'  values  or  large  molecular  weight  discrepancies  require  higher  values  of  X 
Accurate  determination  of  X  values  usually  require  a  full  boundary  layer  solution 
study . 
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SECTION  3 


USER’S  GUIDE 
ASTHMA  PROGRAM 

This  section  of  the  report  summarizes  the  input  requirements  of  the 
ASTHMA3  program  in  a  form  convenient  for  users  of  the  program  and  provides  spe¬ 
cific  instructions  for  program  operation. 

3.1  DESCRIPTION  OF  INPUT 

The  input  to  the  Axi-Syrometric  Transient  Heating  and  Material  Ablation 
Program  can  conveniently  be  considered  as  having  nine  parts,  six  of  »?hich  must 
be  present  for  each  run;  the  remaining  three  are  optional,  depending  on  various 
input/output  choices  and  on  the  heating  option  chosen.  The  individual  nine  parts 
are  described  in  the  nine  following  sub-sections. 

3.1.1  Title  and  Heading  Information 

The  first  three  cards  of  the  data  deck  are  used  to  transmit  title  and 
heading  information  to  the  output.  The  first  72  columns  of  each  of  these  cards 
may  be  used  for  the  title,  the  alphamei'ic  information  in  columns  61  through  77 
of  the  third  card  being  used  as  a  page  heading  on  all  pages  after  the  first. 

3.1.2  General  Problem  Constants 

One  card  provides  the  program  with  certain  general  constants. 


Column 

Format 

Data 

Units 

1-3 

13 

Largest  m  index  in  nodal  network,  number  of 
rows  of  boxes  {not  corners) 

— 

4-6 

13 

Largest  n  index  in  nodal  network,  number  of 
columns  of  boxes  (not  corners) 

— 

7-12 

E6.4 

Initial  value  of  problem  time 

sec 

13-18 

E6.4 

Final  value  of  problem  time 

sec 

19-24 

E6.4 

Output  time  interval  if  greater  than  zero. 

If  blank,  the  variable  output  interval 
option  is  used  (see  Section  3.1.3). 

sec 
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Column  Format  Data  Units 

71  II  0  -  Back  shifted  nodes  - 

1  -  Centered  nodes  - 

72  II  Property  look-ups  and  formulation  of  con¬ 

ductances  and  capacities, 

0  -  Uses  old  AATT  program  procedures  - 

1  -  Uses  more  efficient  procedures  written  - 

specifically  for  ASTHMA3  {not 
operational  in  some  editions) 

73  II  0  -  Thermal  conductance  calculations  based 

on  rectangular  nodal  boxes 

1  -  Thermal  conductance  calculations  using 
sine  correction;  improves  accuracy  for 
nonorthogonal  nodal  boxes 


3.1.3  Special  Output  Interval  Specification  (Optional) 


A  card  with  the  following  information  is  required  if  columns  19-24  of 
the  card  described  in  Section  3.1.2  are  blank. 


Column 

Format 

Data 

Units 

1-10 

F10. 0 

First  output  time  interval,  applied  at 
initial  time 

sec 

11-20 

F10.0 

Transition  time  between  first  output  time 
interval  and  second  output  time  interval 

sec 

21-30 

F10.0 

Second  output  time  interval 

sec 

31-40 

F10.0 

Transition  time  between  second  output  time 
interval  and  third  output  time  interval 

sec 

41-50 

F10.0 

Third  output  time  interval 

sec 

51-60 

F10. 0 

Transition  time  between  third  output  time 
interval  and  fourth  output  time  interval 

sec 

61-70 

F10.0 

Fourth  output  time  interval 

sec 

71-80 

F10.0 

Not  used 

3.1.4  Special  Punched  Temperature  Output  Specification  (Optional) 


If  KSTRP  on  the  card  described  in  Section  3.1.2  (column  63)  is  2,  up  to 
eight  special  output  times  are  to  be  provided  for  the  punched  output  on  a  single 
card.  See  output  description  for  further  information. 


Preceding  page  blank 
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Column 


Format 


Data 


Units 


1-10  8F10.0  Up  to  eight  output  times  for  punched  out-  sec 

put  data 

11-20 

21-30... 

3.1.5  Kodal  Net  Corner  Coordinate  Cards 

This  set  of  cards  provides  the  nodal  box  comer  coordinates  and  certain 
nodal  center  locations  to  the  program.  For  m  rows  and  n  columns,  there  will  be 
(ra+1)  (n+1)  corner  coordinate  cards.  It  is  advisable  to  run  the  corner  coordi¬ 
nate  cards  through  any  convenient  plot  program  to  check  the  nodal  layout  before 
any  heat  transfer  calculations  are  made. 

Corner  cards  must  be  entered  in  order  beginning  at  m=l,  n=l  and  continu¬ 
ing  up  the  first  column  of  corners  (n=l)  for  all  m's  until  column  is  completed, 
returning  to  bottom  of  next  column  m=l,  n=2,  and  so  on,  as  indicated  in  Fig.  17. 


Heated  Surface 


Figure  17.  Sketch  of  Nodal  Corner  Input  Order 

Corners  m,  n  representing  boxes  m,  n  which  have  no  material  (null  nodes 
entered  to  complete  the  mesh  lattice)  may  be  blank  entries,  but  a  card  must  al¬ 
ways  be  included.  There  must  be  (m+1) (n+1)  corder  cards. 

Vresent  dimensions  allow  300  nodal  boxes,  with  a  maximum  of  40  columns. 


Column  Format 


Data 


Units 


1-7  E7.5  Corner  radius,  rQ  inches 

8-1 4  E7.5  Corner  axial  location,  inches 

*The  user  should  keep  in  mind  the  discussion  of  Section  2.2. 2.1.  The  m,n 
network  may  be  located  completely  independently  of  the  physical  r-Z  coordinate 
scheme.  The  heated  surface  is  at  the  "top"  of  a  column,  but  "top"  merely 
means  the  last  coordinate  card  read  in  a  given  column. 
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3.1.6  Nodal  Data  Cards 


This  set  of  cards  identifies  the  materials  in  the  nodal  boxes,  flags  the 
time  tables  and  surface  thermochemistry  tables  to  be  associated  with  the  nodes, 
and  specifies  radiation  view  factors.  For  an  m  by  n  mesh  network  mn  cards  must 
appear,  blank  cards  being  entered  for  null  (no  material)  nodes.  Cards  begin  at 
1,1  and  proceed  to  m,l,  followed  by  1,2  to  m,2,  and  so  on. 


Column 


Format 


Data 


II  Material  number  of  material  .in  this  box, 

refers  to  numbered  material  properties 
tables.  Blank  for  null  node 

II  Enter  1  if  used  for  stability  limit  cal¬ 

culation  of  time  step.  Leave  blank  if 
this  node  can  safely  be  emitted  from 
stability  considerations  (economy  measure) 

II  Enter  1  if  this  node  is  at  the  heated 

surface,  enter  2,  3,  or  4  to  that  side 
assigned  to  backwall  boundary  conditions 
(see  Figure  18  below) ,  otherwise  blank; 
backwall  assignment  rules  are  discussed 
below 
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11  Enter  number  of  surface  thermochemistry 
table  which  is  to  be  used  for  energy 
balance  calculations  on  this  node  if  it 
is  or  becomes  a  surface  node  (see  the 
discussion  immediately  below  and  Sec¬ 
tion  3. 1.8. 3  below) 

12  Enter  number  of  time-function  table 
which  is  to  be  used  for  boundary  condi¬ 
tions  for  this  node  if  it  is  or  becomes 
a  surface  node  or  if  it  is  a  backwall 
(or  side)  node  which  has  been  flagged 
by  2,  3,  or  4  in  column  3.  If  it  is  a 
flagged  backwall  (or  side)  node,  a  zero 
in  column  6  dictates  that  the  heat 
transfer  coefficient  defined  in  columns 
49-54  of  the  data  card  described  in  Sec¬ 
tion  3.1.2  will  be  used.  Other  rules 
on  backwall  assignments  are  discussed 
below. 


7-12 

E6.4 

Leave  blank  (reserved  for  future  use) 

13-18 

E6.4 

Initial  temperature 

of  this  node 

19-24 

E6.4 

Interface  (contact) 

resistance 

at  top  of  node  (m+1  plane) 


Units 


°R 

ft2sec°R/Btu 
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Column 

Format 

Data 

Units 

25-30 

E6.4 

Interface  (contact)  resistance  at  right 
side  of  node  (n+1  plane) 

ft2sec°R/Btu 

31-36 

E6.4 

Option  1  radiation  view  factor  for  this 
node  if  it  is  or  becomes  a  surface  node 
or  is  a  side  or  back-wall  side 

- - 

37-42 

E6.4 

Option  3  radiation  view  factor  for  this 
node  if  it  is  or  becomes  a  surface  node 

— 

Assignments  and  options  for  nodes  require  some  amplification  over  and 
above  the  summaries  given  in  the  brief  column-by- column  descriptions  given 
above: 

Surface  Thermochemistry  Table  Assignment  (Column  4) 

Surface  thermochemistry  tables  are  described  in  more  detail  in  Section 
3.1.9  below.  Each  table  is  made  for  one  pressure  and  one  ablating  material. 
The  tables  are  distinguished  during  the  reading  process  by  the  pressures.  The 
pressure  is  noted  for  each  card,  A  change  in  pressure  signifies  the  start  of 
a  new  table.  The  pressure  sets  or  tables  are  numbered  as  they  are  read  in  be¬ 
ginning  with  number  one.  It  is  this  number  which  is  referred  to  in  the  column 
4  punch. 

Back  Wall  Boundary  Condition  Assignments 

The  backwall  boundary  condition  may  be  applied  to  any  of  three  surfaces 
of  a  node  (other  than  the  "top  of  the  column"  side)  and  to  more  than  one  node 
in  a  column  (excepting  only  the  node  at  the  heated  surface) .  The  backwall 
boundary  condition  may  vary  from  node  to  node  and  with  time  and  will  allow  an 
assigned  temperature  option.  The  following  notes  describe  the  ground  rules 
for  applying  backwall  boundary  condition : 

a.  A  node  is  exposed  to  the  backwall  boundary  condition  by  a  2,  3,  or 
4  punch  in  the  "side  heated"  slot  of  the  nodal  data  card,  where  the 
sides  are  numbered  as  shown  in  Figure  18 


m  index 
(up  a  column) 


1 


4 


(nodal 

box) 


2 


3 


n  index 
(along  a  row) 


Figure  18.  Definition  of  the  Numbering  System  for  the  Sides  of  a  Nodal  Box 
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b.  Nodes  not  assigned  as  "backwall  nodes"  will  be  insulated  on  sides  2, 

3,  and  4  if  these  do  not  communicate  with  other  nodes  (i.e.,  are 

"exposed") 

c.  Backwall  nodes  may  be  assigned  temperature  (Option  2)  or  be  convec- 

tively  and/or  radiatively  heated  or  cooled  (Option  1) .  See  Section 

2.3  for  discussion  of  Options  1  and  2. 

d.  For  Option  1  backwall  nodes: 

(1)  The  reservoir  temperature  with  which  all  such  nodes  communicate 
is  the  constant  value  given  in  column  61-66  of  the  first  numeri¬ 
cal  data  card  (see  Section  3.1.2  above). 

(2)  The  convective  transfer  coefficient  will  be  obtained  from  the 
time  table  assigned  to  this  node  (in  which  table  the  coeffi¬ 
cient  h  in  Btu/ft2sec°F  is  entered  in  columns  31-40  of  each 
card  otherwise  used  for  p  ueCH) ;  however,  if  no  time  table  is 
assigned  to  a  backwall  node,  h  will  be  obtained  from  the  con¬ 
stant  value  entered  in  columns  49-54  of  the  first  numerical 
data  card;  zero  values  of  h  in  an  Option  1  time  table  are  not 
allowed;  radiation  heat  flux  and  recovery  enthalpy  values  in 
these  time  tables  are  ignored. 

(3)  The  radiation  view  factor  will  be  taken  as  the  Option  1  view 
factor  assigned  to  the  node. 

(4)  The  emissivity  or  emittances  of  the  exposed  backwall  surface 
will  be  taken  from  the  material  properties  table  values  of  e 
for  the  material  assigned  to  this  node;  however,  if  this  e  is 
zero,  e  will  be  obtained  from  the  constant  value  entered  in 
columns  55-60  of  the  first  numerical  data  card  (see  Section 
3.1.2  above) . 

(5)  Backwall  convective  and  radiative  transfer  reduces  the  stabil¬ 
ity  time  step  for  a  node;  high  values  of  h  should  be  avoided; 
see  Section  2.2.4  above. 

e.  For  Option  2  backwall  nodes: 

(1)  The  assigned  temperature  will  be  obtained  from  the  time  table 
assigned  to  this  node;  if  no  time  table  is  assigned,  the  node 
cannot  be  treated  under  Option  2. 

(2)  No  recession  may  be  input  in  time  tables  for  backwall  nodes, 
such  nodes  never  recv.^e. 

f.  Backwall  nodes  may  never  switch  heating  options. 
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g.  If  erosion  exposes  a  back  or  sidewall  node  to  the  heated  surface 
boundary  condition,  no  provision  has  been  made  for  an  automatic 
smooth  reinterpretation.  Entries  in  the  time  table  will  be  inter¬ 
preted  in  the  non-backwall  sense  (see  Section  3.1.8  below) .  In 
Option  1,  this  will  lead  to  peculiar  results. 

3.1.7  Materials  Properties  Tables 

One  table  of  temperature  dependent  material  properties  is  input  for  each 
material  appearing  in  the  nodal  network.  A  lead  card  gives  the  total  number  of 
materials  used.  The  tables  are  not  numbered  and  are  presumed  to  be  encountered 
in  ascending  numerical  order,  beginning  at  one  and  including  no  gaps  or  omitted 
table  numbers.  The  main  ablating  material  must  have  material  number  and  table 
number  one.  There  must  be  at  least  two  temperature  entries  in  each  table. 
Present  dimensions  allow  six  material  property  tables,  with  15  temperature  en¬ 
tires  in  each  table. 

First  Card 


Column 

Format 

Data 

Units 

1-2 

12 

Total  number  of  materials 

— 

Tabular  Entry  Cards 

Column 

Format 

Data 

Units 

1-2 

12 

Flag.  Last  entry  in  table  has  -1;  blank 
for  other  cards 

— 

3-8 

E6.4 

Temperature 

°R 

9-14 

E6.4 

Density.  Enter  in  first  card  only;  will 
not  be  a  function  of  temperature 

lb/ft3 

15-20 

E6.4 

Specific  heat  c^ 

Btu/lb°R 

21-26 

E6.4 

Thermal  conductivity  along  rows  (n-direc- 
tion,  fixed  m) 

Btu/ft  sec°R 

27-32 

E6.4 

Emissivity  (emittance) 

— 

33-38 

E6.4 

Thermal  conductivity  along  columns 
(m-direction ,  fixed  n) ,  if  blank  will  be 
taken  as  equal  to  n-direction  conductivity 

Btu/ft  sec°R 

3.1.8 

Functions-of- 

-Time  Tables  ("Heating"  Tables) 

A  series  of  tables  provides  heated  surface  boundary  condition  informa¬ 
tion  as  functions  of  time  and  in  some  cases  back  or  sidewall  boundary  conditions 
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as  well-  The  input  nodes  for  the  two  cases  are  slightly  different,  but  in 
general  the  same  input  format  is  followed. 

The  various  time  tables  are  independent,  and  each  nodal  box  in  the  in- 
depth  layout  has  been  given  the  number  of  a  time-table  for  use  if  and  when 
that  node  becomes  a  surface  node,  or  {in  some  cases) it  is  a  back -wall  node. 

Three  general  types  of  boundary  condition  sets  are  available,  referred 
to  here  as  Options  1,  2,  and  3.  Section  2.3  above  provides  further  descrip¬ 
tion. 

The  three  options  are: 

Option  1  -  General  convection-radiation  heating  with  coupled  mass 
transfer,  including  the  effects  of  unequal  heat  and 
mass  transfer  coefficients  (non-unity  Lewis  number)  and 
unequal  mass  diffusion  coefficients. 

Option  2  -  Specified  surface  temperature  and  surface  recession  rate 

(Surface  recession  rate  is  not  allowed  for  backwall  nodes). 

Option  3  -  Specified  radiation  view  factor  and  incident  radiation 
flux,  as  functions  of  time,  for  a  stationary  surface 
("cooldown"  option).  This  specific  option  is  not  avail¬ 
able  for  backwall  nodes,  but  the  same  physics  can  be 
modeled  with  an  Option  1  table. 

Options  1  and  3  are  surface  energy  balance  options;  Option  2  does  not,  of 
course,  use  an  energy  balance. 

For  heated  surface  boundary  conditions,  each  time  table  may  have  a  se¬ 
quence  of  heating  options;  thus  it  is  convenient  to  think  of  each  time  table 
as  consisting  of  a  number  of  sub-tables,  each  sub-table  representing  one  option. 
The  switch  from  one  option  to  another  requires  a  repeated  time  entry,  the  first 
card  representing  the  last  entry  of  the  earlier  table  and  the  second  card  rep¬ 
resenting  the  start  of  the  next  table. 

The  total  number  of  time  values  in  each  time-table  is  limited  to  25. 

Each  sub-table  must  have  at  least  2  entries,  hence  the  number  of  sub-tables 
in  a  table  cannot  exceed  12.  (The  most  common  problem  has  only  2,  represent¬ 
ing  an  Option  1  or  Option  2  calculation  followed  by  cooldown,  Option  3.)  Time 
entries  must  be  in  increasing  order. 

Backwall  time  tables  may  not  have  option  switches. 

The  tables  are  presumed  to  be  numbered  sequentially  beginning  with  table 
number  1.  There  may  be  ten  different  time  tables.  The  format  for  the  time¬ 
tables  is  as  follows: 
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Column  Format  Data 

1-2  17:  Flag,  nominally  blank,  punched  to  indi¬ 

cate  tl.e  last  card  of  each  time  tnble. 
Punch  +1  for  last  card  of  last  time 
table,  -1  for  last  card  of  any  preceding 
time  tables 

3-10  F8.2  Time  (independent  variable) 

11-20  F10.2  •  Surface  (or  side  1)  heating 

Option  Is  Recovery  enthalpy,  relative  to 
the  same  chemical-base  state  as  used  with 
the  heats  of  formation 

Option  2:  Surface  temperature 

Option  3:  Blank 

•  Backwall  (or  sides  2,  3,  or  4)  heating 
Option  1:  Blank 

Option  2:  Backwall  node  temperature 


Pnits 

sec 

Btu/lb 

°R 

°R 


21-30  F10 . 5  Option  1:  Radiant  energy  flux  to  the 

surface  Btu/ft^  sec 

Option  2:  Surface-recession  measured  along 
nodal  column  centerlines  (must  be  zero  for 
backwall  application)  mils 

Option  3:  Radiant  energy  flux  to  tn®  sur-  ,  , 

face  Btu/ft2sec 


31-40  F10.5  •  Surface  (or  side  1)  heating 

Option  1:  heat  transfer  coefficient  lb/ft2sec°F 

Option  2 :  Blank 

Option  3:  Must  be  blank 

•  Backwall  (or  sides  2,  3,  or  4)  heating 

Option  1:  Heat  transfer  coefficient  Btu/ft  sec°F 

Option  2 :  Blank 

41-50  F10.5  Option  1:  Pressure  (not  used,  information 

only)  atm 

Option  2:  ad  lib 

Option  3:  ad  lib 

51-60  F10.5  Blowing  rate  parameter  A  (see  Section 

3.1.2) 


61-70  F10.5  Heat  rate  q  .  into  a  surface  or  backwall  node 

(Options  1  and  3;  see  Figure  15  and  Equations 

(34),  (35),  (42))  Btu/sec 
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3.1.9  Surface  Thermochemistry  Data  (Option  1  Only) 

3 . 1 . 9 . 1  Introduction 

Problems  involving  use  of  the  surface  thermochemistry  option  (Option  1) 
at  the  heated  surface  require  the  input  of  an  array  of  surface  mass  and  energy 
data  particular  to  the  option  and  the  material  being  analyzed.  These  data  in¬ 
clude  the  specification  of  the  ratio  of  mass  transfer  coefficient  to  heat  trans¬ 
fer  coefficient,  the  radiation  view  factor,  and  surface  thermochemical  data. 

The  cards  containing  this  input  are  described  below.  Problems  not  involving 
Option  1  calculations  do  not  need  the  surface  equilibrium  data  deck.  The  pro¬ 
gram  can  be  instructed  to  use  the  surface  tables  from  the  preceding  problem 
(if  any) ,  in  which  case  additional  surface  tables  are  not  provided.  The  "sur¬ 
face  tables"  then  consist  of  only  the  "lead  card"  described  below.  For  Option 
1  problems  involving  multiple  ablating  materials  see  Appendix  A  for  minor  modi¬ 
fications  to  the  surface  thermochemistry  tables  as  described  here. 


3. 1.9. 2  Ratio  of  Mass  to  Heat  Transfer  Coefficient  and  Kinetics 
Parameter 

A  single  card  serves  to  specify  the  ratio  of  the  convective  mass  trans¬ 
fer  coefficient  to  the  convective  heat  transfer  coefficient  (CM/CH) • 


Column 

Format 

Data 

Units 

1-10 

F10.0 

CM/CH 

— 

11-20 

F10.0 

Not  used 

— 

21-39 

19X 

Blank 

40 

11 

One  punch  will  cause  surface  tables  of  pre¬ 
ceding  job  to  be  used  for  this  job.  Read 
no  more  input.  Otherwise  leave  blank. 

41-49 

9X 

Blank 

— 

50 

11 

Format  flag  (see  Section  3. 1.9. 3. 3). 

0  or  blank  -  standard  ACE  and  GASKET  format 

3  -  CMA  format 

— 

51-56 

11 

Material  number  for  which  cp(T)  will  be 
taken  for  making  surface  thermochemistry 
calculations  using  first  pressure  table 
below  (blank  implies  one) ;  see  Appendix  A 
for  discussion 

F5.0 

Heat  of  formation  to  be  used  in  perform¬ 
ing  surface  thermochemistry  calculations 
using  first  pressure  table  below  (blank 
calls  for  value  entered  in  columns  37-42 
of  first  data  card) ;  see  Appendix  A  for 
discussion 

Btu/lb  at 
536°R 
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Column 


Form*  it 


Data. 


Halt 


57-62 

63-68 

69-74 

75-80 


Same  as  for  51-56  for  second,  third, 
fourth  and  fifth  pressure  tables;  see 
Appendix  A  for  discussion;  leave  blank 
for  simple  ablating  material  problems 


3. 1.9 -3  Surface  Thermochemistry  Table 
3. 1.9. 3.1  Introduction 

This  table  supplies  the  necessary  input  data  for  the  surface  energy 
balance  computations  in  Option  1.  (This  energy  balance  is  discussed  in  Sec¬ 
tion  2.3.) 

Host  commonly  the  deck  of  ca_ds  which  n  ke  up  the  surface  thermochemis¬ 
try  table  is  generated  by  the  Aerotherm  Chemical  Equilibrium  Program  (ACE  or 
GASKET).  The  user's  manuals  for  the  ACE  and  GASKET  programs  describe  this 
table  in  complete  detail  (Refs.  4  and  5) .  On  occasion  the  user  may  desire 
to  construct  his  own  surface  thermochemistry  table,  and  so  the  following  sec¬ 
tions  include  brief  descriptions  of  the  organization  and  forma;.  of  these  tables. 
The  main  emphasis  however  will  be  on  the  coirmun ic at ion  between  the  ACE  or  GASKET 
programs  and  the  ASTHMA  program,  since  this  is  of  the  most  general  interest. 


3 . 1 . 9 . 3 . 2  Edge  Enthalpy  Data 

Equations  (34)  and  (35)  of  Section  2.3  indicate  that  if  diffusion  coef¬ 
ficients  are  not  equal  or  if  the  ratio  Cj-/CH  is  not  unity,  then  the  surface 
energy  balance  requires  data  about  the  edge  gases  of  the  boundary  layer.  These 
data  are  provided  in  special  "edge  tables”  which  precede  each  pressure  section 
of  the  surface  tables  (the  various  sections  of  the  surface  tables  are  described 
in  Section  3. 1.9. 3. 3  below).  The  independent  variables  for  an  edge  table  are 
pressure  and  temperature.  Dependent  variables  are  h  and  the  sum^zf  h-  W  . 

GW  Z— /  16  X 

If  the  format  flag  on  the  preceding  card  is  0,  1,  or  2,  the  edge  enthalpy 
data  are  entered  on  the  cards  as  follows  ("ACE  format") : 


Column 

Format 

Data 

Ui.its 

1-6 

E6.4 

Pressure 

atm 

7-26 

2  OX 

Blank 

— 

27-30 

F4.2 

Unequal  diffusion 

exponent*  y 

— 

This  quantity  is  discussed  in  the  user's  manual  for  the  ACE  program  {Ref.  4) 
and  also  in  References  6-10, 
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Colt 3E3 


Foma" 


Bata 


Units 


32-37 

27.5 

VV=»rr-ti=r^»  f-rnr-fa 

®K 

f®3  if  negative, 
in  which  case 
enthalpies  be¬ 
low  are  Btu/lb) 

38-43 

£6.4 

?Iaused 

44-51 

26. 5 

Surmation  Y* St  h.w 

Z—t  le  a 

cal/gr 

(Btu/lb  if  tem¬ 
perature  is  en¬ 
tered  with  minus 
sign) 

52-59 

28.5 

h 

ev 

cal/gr 

(Btn/lb  if  tem¬ 
perature  i  s  en¬ 
tered  with  minus 
sign) 

60-65 

£6 

Unused 

— 

66 

11 

0  {flag  signifying  that  this  card  is 
part  of  the  edge  gas  table) 

— 

67-78 

2A6 

Problem  identification  (not  read) 

— 

79-80 

12 

Page  number  of  ACS  or  GASSET  output 
listing  containing  the  data  punched 
on  this  card  (not  read) 

” 

For  format 

flag  of  3,  the  following  format  applies: 

Column 

Format 

Data 

Units 

1-8 

F8.5 

Pressure 

atm 

9-16 

J-8.5 

Blank 

— 

17-24 

8X 

Blank 

— 

25-33 

F9.4 

Temperature 

(°R  if  negative 
in  which  case 
enthalpies  be¬ 
low  are  Btu/lb) 

34-38 

F5.3 

Unequal  diffusion  exponent  y 

— 

39-47 

F9.3 

T  , 

Summation  y 

cal/gr 

(Btu/lb  if  tem¬ 
perature  is  en¬ 
tered  with  minus 
sign) 
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Column 

Format 

Data 

Units 

48-55 

F9.3 

Enthalpy  of  edoe  oases  h 

cal/gr 

(Bfcu/lb  if  tem¬ 
perature  is  en¬ 
tered  with  minus 
sicn) 

57-58 

12 

— 1  (flag  signifying  that  this  card 
is  part  of  the  edge  gas  table) 

— 

59-60 

2X 

Blank 

— 

61-66 

A6 

Unused 

— 

67-78 

2P.6 

Problem  identification  (not  read) 

— 

79-80 

12 

Pace  number  on  ACS  or  GASKET  output 
listing  containing  the  data  punched 
on  this  card  (not  read) 

— 

2«ote  that  although  ACE  anc  GASKET  wi.ll  provide  data  decks  using  °K  and 
cal/gr,  in  those  rare  cases  in  which  a  user  wishes  to  supply  his  own  deck  and 
prefers  to  work  in  °R  and  3tu/lb,  he  cay  do  so  simply  by  introducing  a  minus 
sign  as  a  flag  in  front  of  the  temperature  entries. 

The  table  length  is  limited  to  5  pressure  sets  (it  nay  have  only  1  pres¬ 
sure  set)  with  not  core  th^n  25  nor  less  than  3  temperature  entries  in  each  set. 
The  series  of  temperature  values  may  be  different  for  the  edge  table  at  each 
pressure  set.  The  table  is  organized  as  a  series  of  sections,  each  represent¬ 
ing  one  pressure  and  each  preceding  the  corresponding  pressure  group  of  the 
surface  thermochemistry  deck  as  described  below.  The  temperature  entries  with¬ 
in  each  section  must  be  ordered,  either  ascending  or  descending.  Similarly, 
the  pressures  must  be  ordered  either  ascending  or  descending.  (Decks  generated 
by  the  ACE  and  GASKET  programs  will  have  been  automatically  ordered  properly.) 

3. 1.9. 3. 3  Surface  Thermochemistry  Tables 

3. 1.9. 3. 3.1  Description  of  Surface  Thermochemical  Tables 

This  table  comprises  a  series  of  sections.  Each  section  represents  one 
pressure*  and  one  transfer  coefficient  value.  Each  section  consists  of  two 
subsections.  The  second  represents  surface  temperatures  too  low  for  ablation; 


The  word  "pressure"  is  used  as  a  handy  identification  parameter.  In  general, 
of  course,  there  are  two  independent  edge  state  variables,  which  are  properly 
accounted  for  by  the  ACE  program  in  its  chemistry  solution,  as  well  as  a  speci¬ 
fication  of  the  ablating  material.  The  ASTHMA  program  needs  only  to  be  given 
some  "flag"  to  distinguish  the  tables.  Pressure  is  chosen  for  this  purpose, 
and  hence  the  word  "pressure"  is  used  here  as  shorthand. 


in  this  subsection  surface  temperature  is  the  independent  variable-  The  first 
subsection  contains  the  ablating  cases;  here  the  ablation  rate  is  the  third 
independent  variable  and  the  surface  temperature  is  a  dependent  variable - 

Thus  one  table  has  three  independent  variables:  pressure,  transfer  coef¬ 
ficient,  and  either  surface  temperature  or  ablation  rate,  depending  on  whether 
the  surface  temperature  is  high  enough  for  ablation. 

The  table  has  either  two  or  three  dependent  variables,  according  to 
whether  the  surface  temperature  is  high  enough  for  ablation.  Two  dependent 

EX  Ay 

Z-wh.  and  h^,  the  enthalpy  of 

the  wall  gases.  The  third  dependent  variable  is  the  surface  temperature,  but 
it  is  dependent  only  in  those  cases  for  which  the  surface  temperature  is  high 
enough  for  ablation.  (Otherwise,  the  surface  temperature  functions  as  an  inde¬ 
pendent  variable.) 

The  ACE  and  GASKET  programs  generate  separate  groups  for  each  pressure, 
one  at  a  time.  All  these  groups  together  make  up  the  surface  thermochemistry 
deck.  Within  each  pressure  group  the  transfer  coefficient  values  will  be 
ordered.  Within  each  transfer  coefficient  section,  non-zero  ablation  rate  en¬ 
tries  will  be  grouped  ahead  of  the  zero  ablation  rate  entries.  The  non-zero 
ablation  rate  entries  will  not  be  ordered  in  any  particular  way  on  the  abla¬ 
tion  rates;  any  necessary  ordering  is  made  automatically  by  the  ASTHMA  program 
as  it  reads  in  the  data.  The  zero  ablation  rate  entries  are  ordered  with  de¬ 
scending  temperatures. 

(Users  providing  their  own  thermochemistry  decks  must  ensure  that  the 
transfer  coefficients  are  ordered,  but  the  ordering  may  be  either  ascending  or 
descending  in  each  case.  Within  each  transfer  coefficient  section,  the  non¬ 
zero  ablation  rate  entries  and  the  zero  ablation  rate  entries  may  be  shuffled 
together  but  the  zero  ablation  rate  entries  must  be  encountered  in  order  of 
descending  temperature  during  the  reading  process.  Non-zero  (ablating)  abla¬ 
tion  entries  need  not  be  ordered.  The  ablation  rate  for  non-ablating  (zero 
ablation  rate)  entries  need  not  actually  be  zero,  but  may  not  exceed  an  abla¬ 
tion  rate  in  the  ablation  section.  These  cards  are  identified  as  zero  abla¬ 
tion  rate  cards  by  a  unity  flag  in  column  66,  as  described  in  the  format  speci¬ 
fication  below.  (The  flag  is  a  zero  in  columns  57-58  for  the  CMA  format.) 

The  number  of  pressure  groups  may  not  exceed  5  (and  may  be  only  1) ;  the 
number  of  transfer  coefficient  values  in  each  pressure  group  may  not  exceed  8 
but  may  be  only  one.  The  sequence  of  transfer  coefficient  values  need 
not  be  the  same  in  the  different  pressure  sections.  Within  each  transfer  coef¬ 
ficient  section  the  number  of  ablation  rate  entries,  including  the  zero 
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ablation  rate  (independent  surface  temperature)  cards,  nay  not  exceed  24  and 
nay  not  be  less  than  2.  The  series  of  ablation  values  3’  nay  be  unique  for 
each  section. 

The  °R-Btu/lb  option  described  for  the  edge  tables  in  Section  3. 9. 3. 2 
nay  be  used  for  these  tables  also. 

3. 1.9. 3. 3. 2  Standard  ACE  Format  (0  Punch) 

If  the  format  flag  given  in  column  50  of  the  card  described  in  Section 

3. 1.9. 2  is  0,  the  card  format  for  the  surface  equilibrium  data  as  produced  by 
GASKET  is  as  follows: 


Column 

Format 

Data 

Units 

1-6 

E6.4 

Pressure 

atm 

7-12 

E6.4 

Gas  rate  ra  /p  u  C..  =  B*  (not  used  by 
ASTHMA3)  9  e  9 

— 

13-20 

E8.5 

Ablation  rate  mc/PeueCH  =  b'c 

— 

21-26 

E6.4 

Transfer  coefficient  peueCH 

lb/ft 2 sec 

27-30 

F4.2 

Unequal  diffusion  exponent  y 

— 

31-37 

E7.5 

Surface  temperature 

°K 

(°R  if  negative 
in  which  case 
enthalpies  be¬ 
low  are  Btu/lb) 

38-43 

E6.4 

Bl  =  m  /p_u  C  (not  used  in  ASTHMA3) 

I  r  £  G  6  H 

— 

44-51 

E8.5 

T 

— -  W 

Summation  >  'Z7  h. 

Z-/  iw  l 

cal/gr 

(Btu/lb  if  tem¬ 
perature  is  en¬ 
tered  with  minus 
sign) 

52-59 

E8.5 

Enthalpy  of  wall  gases  hw 

cal/gr 

(Btu/lb  if  tem¬ 
perature  is  en¬ 
tered  with  minus 
sign) 

60-65 

A6 

Chemical  symbol  of  surface  species. 

(ACE  and  GASKET  print  such  symbols 
arranged  alphabetically  and  trun¬ 
cated  from  right  end  if  necessary) 

66 

11 

1  for  assigned-temperature  entries  in 
the  equilibrium  program  (no  ablation) ; 

>  1  for  surface  thermochemistry  with 
ablation  (temperature  is  dependent) 

Z - 

Provided  only  by  GASKET,  not  by  ACE. 
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Column 


Format 


Data 


Onits 


67-78  2A6  Problen  identification  (not  read) 

79-86  12  Page  number  on  ACE  output  listing  contain¬ 

ing  the  data  punched  on  this  car-'  (not 
read) 


3. 1-9. 3. 3*3  CMA  Format  (3  Punch) 

A  format  flag  of  3  calls  the  CMA  format-  This  is  an  optional  ACE  and 
GASKET  output  format  intended  for  use  with  th1:  CMA  program  (Ref-  10)  .  Since  so 
many  graphite-in-air  surface  tables  already  :xist  with  this  format  it  is  useful 
to  have  it  as  an  ASTHMA3  input  format  also-  In  this  format  the  transfer  coef¬ 
ficient  is  assumed  by  ASTHMA3  to  be  in  whcfc  is  the  slot  of  the  CMA  input. 

An  ACE  problem  with  kinetics  output  in  tie  CMA  card  format  will  lose  the  ki¬ 
netics  information,  however,  since  therr  is  no  provision  for  it  in  the  ACE 
output  statement-  Hence,  for  ACE  outp*.;  this  format  will  probably  not  be  used 
for  problems  with  kinetic  control,  an'-  the  transfer  coefficient  values  entered 
in  the  slot  will  be  fake  numbers-  Output  from  GASKET  will  have  the  transfer 
coefficient  punched  properly. 


Column 

Format 

Data 

Units 

1-8 

F8.5 

Pressure 

atm 

9-16 

F8.5 

Transfer  coefficient* 

lb/ft 2 sec 

17-24 

F8.5 

Char  rate  mc/Peuecn 

— 

25-33 

F9.4 

Surface  temperature 

°K 

( °R  if  negative 
in  which  case 
enthalpies  be¬ 
low  are  Btu/lb) 

34-38 

F5.3 

Unequal  diffusion  exponent  y 

— 

39-47 

F9.3 

T 

Summation  Y^Z*  h. 

jL-J  1W  1 

cal/gr 

(Btu/lb  if  tem¬ 
perature  is  en¬ 
tered  with  minus 
sign) 

48-56 

F9.3 

Enthalpy  of  wall  gases  hw 

cal/gr 

(Btu/lb  if  tern- 

perature  is  en¬ 
tered  with  minus 
sign) 


-k 

Not  provided  by  ACE;,  provided  properly  by  GASKET. 
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Colinn 


Units 


Format 


Data 


57-58  12  0  for  assigned-temoerature  entries  in  the 

thermochemistry  table  (no  ablation) ;  >  0 
for  surface  thermochemistry  with  ablation 
(temperature  is  dependent) 

59-60  2X  Blank 


61-66  A6  Chemical  symbol  of  surface  species.  (ACE 

and  GASKET  programs  print  such  symbols 
arranged  alphabetically  and  truncated  from 
right  end  if  necessary.) 

67-78  2A6  Problem  identification  (not  read) 


79-80  12  Page  number  on  ACE  or  GASKET  output  list¬ 

ing  containing  the  data  punched  on  this 
card  (not  read) 


3. 1.9. 3. 4  Termination  Card 

The  surface  equilibrium  data  deck  must  be  terminated  by  a  single  blank 
card.  Outpuc  decks  of  the  ACE  and  GASKET  programs  may  not  have  such  a  card,  in 
which  case  the  user  must  supply  it. 


3. 1.9. 3. 5  Assembled  Thermochemical  Peck 

Figure  19  shows  a  picture  of  an  assembled  thermochemical  data  deck  for 
several  pressures. 


3.2  PROGRAM  OUTPUT 

3.2.1  Input  Data  (Except  Surface  Thermochemical  Tables) 

Program  output  begins  with  an  output  of  the  input  title  and  heading 
information,  output  interval  specifications  and  general  program  constants, 
nodal  data  (including  any  computed  nodal  center  locations) ,  material  property 
tables,  and  time  dependent  boundary  conditions  table. 

All  this  output  is  fully  labeled  and  is  printed  exactly  as  input  by  the 

user. 


3.2.2  Surface  Thermochemistry  Tables 
3. 2. 2.1  Edge  Enthalpy  Table 

If  there  is  an  edge  enthalpy  table,  it  is  output  exactly  as  input. 
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Figure  19.  Sketch  of  Surface  Thermochemistry  Table  Make-Up 


3. 2. 2. 2  Surface  Thermochemistry  Table 


This  table  is  output  re-ordered  with  increasing  ablation  rates  in  each 
section  and  with  a  new  computed  term  in  place  of  the  input  enthalpy  terms.  For 
each  entry  in  the  surface  thermochemistry  tables  the  program  computes  the 
quantity 


T 

Z*.  ) h.w 
1W  l 


B'h  + 
c  w 


m  h 
c  c 


peueCM 


qchem  _  Chem 
peueCM  prod 


(42) 


and  then  output  it  as  the  dependent  variable  of  interest  in  the  output  surface 
equilibrium  table.*  This  term  occurs  directly  in  the  surface  energy  balance 
(Equations  (34)  and  (35)  of  Section  2.3.3  above)  and  is  a  useful  diagnostic 
quantity.  In  the  output,  this  quantity  is  labeled  CHEM  PROD.  It  has  the  units 
(Btu/ft2sec)/(lb/ft2sec)  =  Bt.u/lb  and  can  be  loosely  thought  of  as  the  chemical 
energy  release  per  pound  of  PeueCy,  where  PeueCM  can  be  thought  of  as  the 
"scrubbing  flux"  or  "Reynolds  flux." 

In  the  output  table,  the  quantity  PeuecM  is  abbreviated  to  CM. 


3-2.3  Regular  Output 
3. 2. 3.1  Introduction 

At  each  output  interval,  as  specified  by  the  user,  the  program  prints 
out  the  current  values  of  ablation  rates,  nodal  temperatures,  and  other  sup¬ 
plementary  information,  as  described  below. 


3. 2. 3. 2  General  Information 


The  first  line  of  output  gives  the  current  problem  time.  The  second 
line  shows  the  current  values  of  a  number  of  miscellaneous  quantities: 


Heading 
QTOT,  SUR 

QTOT,  INT 

CNSV  ENER 


Total  energy  input  to  the  surface 
SqconddAd6'  Btu 

Total  energy  rise  of  subsurface  mate¬ 
rials,  Btu 

Conservation  of  energy  check,  ratio 
QTOT , SUR/ QTOT , INT 


If  edge  tables  ,are  omitted,  only  the  last  three  terms  on  the  left  are  included. 
T.  a  mechanical  removal  energy  term  has  been  omitted  from  this  equation  since 
it  is  not  presently  considered  by  ASTHMA3 . 
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Heading 


CRNODE 

ITER 

NODE  D-TIME 

ACT  D-TIME 

3. 2. 3. 3  Miscellaneous  Surface  Data 

A  single  line  gives  a  set  of  da 
Each  line  presents  the  following  data: 


m,  n  indices  of  node  yielding  smallest 
stability  limited  time  step 

Total  number  of  computational  cycles 
(steps)  required  to  reach  the  indicated 
value  of  output  time.  Mote  chat  the 
initial  time  is  numbered  as  1. 

Stability  limited  time  step  for  critical 
node ,  sec 

Actual  current  time  step;  may  be  smaller 
than  stability  time  step  due  to  high  sur¬ 
face  recession  or  output  matching 


for  each  point  on  the  heated  surface. 


Heading 


ROW 

Row  number  of  node  below  this  surface 
point  (remember  that  surface  points  are 
distinct  from  nodal  points) 

CCL 

Column  number  of  this  surface  point 

OPTN 

Current  problem  (heating)  option  at  this 
surface  point  (1,  2,  or  3) 

SURF  ITER 

Number  of  iterations  required  for  the 
surface  energy  balance  during  the  pre¬ 
vious  cycle 

SURF  TEMP (R) 

Surface  point  temperature  in  degrees 
Rankine 

H  EDGE 

The  input  local  recovery  enthalpy  as 
determined  by  linear  interpolation  in 
the  f unctions-of-time  table 

H  WALL 

Option  1:  Enthalpy  of  edge  gases  at  the 
wall  temperature,  h 

Option  2 :  Blank 


Option  3:  Blank 

B  PRIME  TOT  Parameter  B'/p^^  =  (B£  +  B£c)/p  ueCM 

(ASTHMA3  requires  B£  =  0) 
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Heading 


HEAT  COEFF 


CH/CHO 


PRESSURE 


RADIUS 

Z 


The  current  value  of  the  convective  heat 
transfer  coefficient 

peueCH 

as  linearly  interpolated  in  the  functions- 
of-time  table  and  corrected  for  the  effect 
of  transpiration  (blowing)  according  to 

C/C„  =  £/(e^-l)  where  £  =  2Xm/C„  and 
H  H0  H0 
C„  is  the  heat  transfer  coefficient  be- 
H0 

fors  being  corrected  for  blowing 

The  ratio  C„/CtI  ,  indicating  the  amount 

H  Ho 

of  blowing  correction 

Current  value  of  pressure  for  this  column 
is  interpolated  in  functions-of-time 
tables 

Current  radial  coordinate  of  this  surface 
point  (inches) 

Current  axial  location  of  this  surface 
point  (inches) 


3. 2. 3. 4  Surface  Kate  Quantities 

A  series  of  lines,  one  for  each  surface  point,  gives  the  current  reces¬ 
sion  rates,  mass  loss  (ablation  rates,  and  surface  energy  flux  rates. 


Heading 

ROW 

COL 

RECESSION  RATES 
CENTER  LINE 

NORMAL 


MASS  RATES 
MDOT  TOTAL 

MDOT  TCHEM 


Row  number  of  node  below  this  surface 
point 

Column  number  of  this  surface  point 

Rate  of  movement  of  surface  point  down 
column  center  line  of  surface  nodal  box 
(mils/sec) 


Rate  of  movement  of  surface  in  direction 
normal  to  tangent  to  surface  at  this 
surface  point  (mils/sec) 


Ate  + 

(ablation) , 


total  local  mass  loss  rate 
lb/ft3  sec* 


m 


tc 


-  « 

In  ASTHMA3,  m  must  be  zero 
rZ 
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Heading 


SURFACE  ENERGY  FLUX  RATES 

This  group  of  terms  all  appear  in  the  surface  energy  balance. 
Equation  (35),  and  have  the  units  Btu/ftssec. 

CONVECTED  IN  qsen  =  PeVVHr  -  h^  ) 

W 

grad  =  awqrad 

grad  =  FewOTw4 
out 

gchem  =  PeueCM(chem  prod) 

CONDUCTION  AWAY  g  , 


RADIATED  IN 
RADIATED  OUT 

CHENICAL  GENERATION 


3. 2. 3. 5  Surface  Time  Integrated  Quantities 

These  lines  correspond  exactly  to  the  line  of  surface  rate  quantities 
described  above.  Here  the  quroitities  are  integrated  over  time  and  the  column 
surface  area,  so  that  the  units  of  these  quantities  are  "Btu  for  this  column  n' 
or  Btu/col,  so  to  speak.  The  area  integration  is  done  since  the  surface  area 
of  a  column  varies  with  time  and  there  is  no  generally  preferred  area  upon 
which  to  base  the  time  integrated  flux  terms. 


3. 2. 3. 6  In-Depth  Data 

This  block  of  data  gives  the  current  temperatures  (degrees  Rankine)  of 
the  subsurface  nodal  points.  The  m,n  index  coordinates  of  the  points  are 
included  for  convenience. 


3. 2. 3. 7  Optional  Punched  Card  Format 

Optional  output  called  for  by  the  KSTRP  flag  of  Section  3.1.2  above  pro¬ 
vides  punched  card  temperature  input  for  various  thermal  stress  codes.  These 
cards  are  punched  for  all  nodal  and  surface  point  temperatures  for  either 
special  punch  times  or  for  all  regular  output  times.  The  card  format  is: 


hew  and  ^ z£ehiW  are  not  available  for  the  output  of  qsen  and  q  h  .  There¬ 
fore  CONVECTED  IN  is  output  as  peueCHHr  (the  cold  wall  heat  flux.  Blowing 
reduction  effect  included)  and  CHEMICAL  GENERATION  is  the  sum  of  all  other 
terms  in  the  surface  energy  balance.  Equation  (34),  except  conduction  and  ra¬ 
diation  quantities.  Since  this  division  of  terms  offers  little  physical  in¬ 
sight,  many  users  prefer  to  use  frozen  edge  tables  even  when  they  are  not 
formally  required. 
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Column 

Format 

Data 

Units 

1-10 

F10.3 

Point  radius  r 

in 

11-20 

F10.3 

Point  axial  coordinate  Z 

in 

21-30 

F10.3 

Point  temperature 

°R 

31-38 

8H 

Units  ININDEGR 

— 

39-41 

13 

Point  row  index  m  (blank  if  a  surface 
point) 

— 

42-45 

IH/,13 

Point  column  index 

— 

46-48 

3H 

Title  MAT  (material) 

— 

49-50 

12 

Material  number  at  this  point 

— 

51-57 

F7.2 

Time 

sec 

58-59 

2H 

Sb,  units  of  time 

— 

60-71 

2A6 

Problem  identification  of  columns  61- 
72  of  third  title  card  (Section  3.1.1) 

— 

72 

IX 

— 

73-75 

13 

Card  number 

— 

76-77 

2H 

OF 

— 

78-80 

13 

Total  number  of  cards  punched  for  this 
time 

— 

3.3 

DUMPS 

3.3.1 

Introduction 

To  prevent  the  execution  of  computations  that  are  wasteful  or  probably 
erroneous,  the  program  provides  certain  emergency  stops  with  dumps  of  diag¬ 
nostic  information. 


3.3.2  Too-Small-Time-Step  Dump 

A  dump  occurs  if  the  program  selects  a  time  step  as  small  as  10” *  sec¬ 
onds.  Computation  ceases  and  the  program  prints  out  the  current  surface  column 
number.  Then  the  program  sets  the  final  problem  time  equal  to  the  current  time, 
and  thus  forces  an  output  of  all  the  current  values  of  the  standard  output 
quantities,  followed  by  a  termination  of  computation. 
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3.3.3  Too-Many-Surface-Iterations  Dump 

The  program  allows  51  iterations  to  find  an  acceptable  surface  energy 
balance  at  any  column  point.  Xf  a  balance  has  not  been  obtained  after  51 
iterations,  the  program  writes  the  diagnostic  message  "Iteration  Stop"  and 
a  block  of  diagnostic  data.  Since  stops  of  this  kind  are  complex  in  nature, 
and  generally  involve  the  surface  thermochemistry  data  deck,  the  dump  data 
must  often  be  communicated  to  the  program  authors  for  analysis. 


Iteration  stops  are  almost  always  due  to  "unfortunate"  variations  in 
the  tabular  energy  quantities  making  up  the  surface  energy  balance.  Hence 
the  iteration  stop  dump  contains  the  51  sets  of  paired  values  of  the  inde¬ 
pendent  variable  (either  temperature  or  in  B^)  and  the  surface  energy  balance 
error  generated  during  the  search  for  a  surface  energy  balance.  The  program 
then  dumps  the  surface  energy  balance  equation  error  (departure  from  zero) 
at  each  tabular  (or  temperature  if  failure  occurred  in  nonablating  section) 
in  the  two  tables  bracketing  the  current  value  of  transfer  ccefficient, 
properly  interpolated  at  the  current  value  of  transfer  coefficient,  for  current 
values  of  CH,  qra(j  ^n»  and  in-depth  conduction  data.  This  information,  plus 
the  iteration  history,  permits  a  rather  complete  pictu.  a  of  the  surface  energy 
balance  error  function  shape,  and  this  in  turn  should  usually  allow  an  accurate 
assessment  of  the  particular  feature  of  the  function  currently  hindering  convergence. 


3.3.4  Unacceptable-Surface-Thermochemistry-Table  Stop 

A  series  of  checks  built  into  the  input  routine  serve  to  detect  common 
errors  in  the  make-up  of  the  surface  thermochemistry  tables.  Discovery  of  an 
error  stops  the  reading  process  and  the  program  prints  o*  '  single  line: 

Bad  Surface  Equilibrium  Table  of  Type  - 

Five  error  types  are  detected  as  follows: 

Type 

0  c  vt  section  of  independent  tempera- 

tuxe  entries  (no  ablation)  is  not  in 
descending  order  in  temperature 

1  Mass  rate  in  no-ablation  subsection 

has  been  set  greater  than  a  mass  rate 
in  the  ablation  subsection.  Although 
the  "zero-mass-rate"  can  differ  from 
zero,  they  must  not  exceed  an  actual 
ablation  rate  entry. 

2  Edge  table  has  been  omitted  even  though 

the  diffusion  coefficients  are  not 
equal  (unequal  diffusion  exponent  £  0) . 
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Type 

3 


4 

c; 


Edge  table  has  been  omitted  even  though 

VCH  *  1 


Inconsistent  unequal  diffusion  exponent 
(this  must  be  uniform  for  all  tables) 

Only  one  or  no  B's  in  a  table 


3.3.5  Sense  Switch  Dumps 

Sense  Switch  1  calls  for  output  every  time  step.  Sense  Switch  2  generates 
extra  diagnostic  dumps  at  every  output,  and  Sense  Switch  3  dumps  some  additional 
surface  table  information  just  prior  to  execution  of  the  main  program  iteration 
loop. 


3 . 4  MISCELLANEOUS 

3.4.1  Running  Time 

On  the  CDC  6600  the  execution  time  for  ASTHMA  is  well  approximated  by 


T(hrs)  =  3  x  10-7  NI  +25%  (41) 

where  N  =  number  of  nodes  and  I  =  number  of  iterations,  provided  that  the 
number  of  iterations  per  output  page  exceeds  about  3  (otherwise  output  time 
begins  to  be  appreciable) .  Iterations  I  may  be  estimated  from 

t  =  Problem  time 

time  step  limit 

and  the  time  step  limit  may  usually  be  estimated  fairly  closely  from  a  small 
number  of  obvious  candidate  nodes: 


time  step  limit  sa  -r— 


for  square  nodes  of  side  i,  and 


J>3 

time  step  limit  4— 


for  "thin"  nodes  of  minimum  side  SL . 

Estimates  for  other  machines  were  made  using  empirical  speed  ratios 
based  upon  similar  Fortran  programs  for  jobs  of  roughly  10  minutes  duration. 
If  we  write  T  =  ANI,  then  we  have 
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Machine  A  Estimated  Error 


(hrs) 

—7 

Philco  2000-212 

5 

X 

10 

±25% 

CSC  33GC> 

5 

X 

io-7 

±50% 

CDC  6600 

3 

X 

10“' 

±25% 

Univac  1108 

3 

X 

10-7 

±50% 

I&M  7094 

12 

X 

10“7 

±50% 

360/40 

50 

X 

io-7 

±50% 

360/44 

25 

X 

io-7 

±50% 

360/50 

25 

X 

io-7 

±50% 

360/65 

5 

X 

io-7 

±50% 

3-4.2  Storage  Requirements 

The  ASTHMA  program  itself  requires  only  about  5200^Q  words  of  storage. 
Associated  subroutines  require  about  1000^Q  words,  flu-t  storage  requirement 
is  for  dimensioned  variables  in  various  commons.  With  present  dimensions, 
these  require  about  15,600-^q  words.  Thus  the  program  is  a  tight  fit  in  a 
32,000  word  machine,  with  present  dimension  values. 

3.4.3  Tape  Requirements 

The  program  uses  nc  scratch  tapes.  It  uses  logical  unit  number  5  for 
input  and  6  for  printed  output,  and  uses  a  punch  tape  through  a  PUNCH  state¬ 
ment. 

3.4.4  Fortran  Deck  Make-Up 

The  ASTHMA  program  is  in  Fortran  IV  and  consists  of  the  following  units 

(1)  ARCAST 

Main  program 

(2)  Subroutine  SURFB 

Surface  energy  balance  procedures 

(3)  Subroutine  LCOUNT 

Counts  lines,  turns  and  numbers  pages 

(4)  Subroutine  LOOK 

Table  look-up  with  linear  interpolation 

(5)  Subroutine  SLOPQ 

Quadratic  curve  fit  and  slope  finder 
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(6)  Subroutine  OGLE 

Table  look-up  with  cubic  curve  fit 

(7)  Subroutine  ORDERD 

Ordering  routine 

(8)  Subroutine  SEQUA 

Orders  according  to  results  of  ORDERD 

(9)  Subroutine  SLOPL 

Linear  curve  fit  and  slope  finder 

(10)  Subroutine  VCOS 

Evaluates  cosine  of  angle  between  surface  normal  and 
nodal  center  line 

(11)  Subroutine  GAP 

Dummy  version  of  routine  which  was  used  to  account 
for  effect  of  radiation  transfer  across  gaps  (no 
longer  used) 
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SECTION  4 


SAMPLE  PROBLEM 

This  section  presents  an  input  listing  and  selected  output  page  listings 
for  a  sample  problem  using  the  ASTHMA3  code. 

The  problem  selected  is  one  of  determining  the  transient  thermal  re¬ 
sponse  of  a  pyrolytic  graphite  coated  substrate  located  at  the  throat  of  a 
rocket  nozzle.  Figure  20  shows  the  nodal  gridwork  used  to  model  the  actual 
configuration . 

An  insulated  boundary  condition  is  imposed  on  the  substrate  backwall. 

The  total  problem  time  is  35.0  seconds  consisting  of  a  29.3  second  firing  and 
a  5.7  second  cooldown.  Output  was  desired  at  every  1  second  interval.  The 
surface  vj.jw  factor  and  the  incident  radiant  energy  flux  for  each  column  dur¬ 
ing  the  firing  (Option  1)  are  equated  to  zero  since  the  propellant  is  trans¬ 
parent  and  the  nozzle  surface  essentially  sees  only  itself.  By  a  similar 
argument  for  cooldown  the  Option  3  radiant  energy  flux  and  the  Option  3  view 
factors  are  set  equal  to  zero.  Figure  20  also  shows  the  locations  and  corre¬ 
sponding  numbers  of  the  surface  thermochemistry/edge  enthalpy  tables  and  the 
functions-of-time  tables.  The  edge  static  pressure  at  the  upstream  thermo¬ 
chemistry/edge  enthalpy  table  (Table  1)  is  56.0  atm  and  37.2  atm  at  the  throat 
(Table  2) . 

The  following  tables  present  a  listing  of  the  problem  input  and  selected 
output  pages. 
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RADIAL  DISTANCE 


AXIAL  DISTANCE,  IN. 


Figure  20.  Sample  Problem  Nodal  Layout 
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OUTPUT,  SAMPLE  PROBLEM  1  (CONTINUED) 
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OUTPUT r  SAMPLE  PROBLEM  1  (CONTINUED) 
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OUTPUT,  SAMPLE  PROBLEM  1  (CONCLUDED) 
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APPENDIX  A 


MULTIPLE  ABLATING  MATERIALS  CAPABILITY 

More  than  one  material  may  appear  on  the  surface  at  a  given  time,  and 
the  materials  exposed  may  change  with  time  (e.g.,  one  material  is  ablated  away, 
exposing  a  new  material) .  Since  the  multiple  ablating  material  feature  is  sel¬ 
dom  used,  the  instructions  for  it  have  been  collected  in  this  appendix. 

Use  of  the  multiple  ablating  materials  capability  in  Option  1  calcula¬ 
tions  requires  special  input  information  for  surface  thermochemistry  table 
assignments.  Each  pressure  table  must  be  assigned  to  a  material  identification 
number,  so  that  the  specific  heat  of  the  proper  ablating  material  can  be  inte¬ 
grated  into  the  surface  thermochemistry.  Furthermore,  a  heat  of  formation  must 
be  input  for  each  ablating  material.  These  assignments  have  all  been  relegated 
to  the  surface  thermochemistry  deck  lead  card  (see  Section  3. 1.9. 2).  Columns 
51-80  of  this  card  are  divided  into  five  fields  each  of  II,  F5.Q  format  which 
are  to  contain  the  following  information: 

II  Material  number  for  which  c^(T)  will  be  taken  for  amalgamation 
into  surface  thermochemistry  calculatons 

F5.0  Heat  of  formation  (Btu/lb  at  536°R)  to  be  used  in  this 
amalgamation 

The  five  material  number  and  heat  of  formation  assignments  will  be  as  assigned 
by  the  ASTHMA3  program  as  needed  to  the  following  pressure  tables  (up  to  five 
in  number).  A  blank  material  number  calls  for  the  use  of  c  (T)  of  material 
number  one  and  of  the  beat  of  formation  entered  in  columns  37-42  of  the  first 
numerical  data  card  (see  Section  3.1.2). 

Since  many  different  nodes  may  be  assigned  to  a  given  thermochemistry 

table  (pressure),  and  any  material  may  appear  in  a  given  node,  then  (since  a 

given  thermochemistry  table  is  assigned  to  a  given  node)  it  is  possible  to  in 

effect  assign  the  same  thermochemistry  table  to  different  materials.  The  ii.put 

conventions  described  above  serve  to  remind  the  user  that  he  must  choose  one 

c  (T)  to  assign,  and  one  heat  of  formation  for  a  given  table.  Therefore,  if 
P 

more  than  one  material  is  assigned  (through  nodal  assignments)  to  a  given 
thermochemistry  table,  these  materials  must,  for  overall  problem  consistency, 
be  thermodynamically  identical  (i.e.,  have  the  same  c^(T)  and  heat  of  formation). 


A-l 


Such  materials  may,  however,  have  different  thermal  conductivities  or  different 
densities . 

Note  that  for  a  problem  with  a  single  ablating  material,  it  is  most  con¬ 
venient  to  number  this  material  as  number  one.  Columns  51-80  of  the  lead  card 
are  then  left  blank,  and  the  heat  of  formation  is  entered  in  columns  37-42  of 
the  first  numerical  data  card. 


